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Abstract 



> 

op Some relevant transport properties of solids do not depend only on the spectrum 

of the electronic Hamiltonian, but on finer properties preserved only by unitary equiv- 
alence, the most striking example being the conductance. When interested in such 



(^ properties, and aiming to a simpler model, it is mandatory to check that the sim- 

^^ pier effective Hamiltonian is approximately unitarily equivalent to the original one, in 

. . the appropriate asymptotic regime. In this paper, we consider the Hamiltonian of an 

^^ electron in a 2-dimensional periodic potential (e.g. generated by the ionic cores of a 

S^ crystalline solid) under the influence of a uniform transverse magnetic field. We prove 

that such Hamiltonian is approximately unitarily equivalent to a Hofstadter-like (resp. 
Harper-like) Hamiltonian, in the limit of weak (resp. strong) magnetic field. The result 
concerning the case of weak magnetic field holds true in any dimension. Finally, in the 
limit of strong uniform magnetic field, we show that an additional periodic magnetic 
potential induces a non-trivial coupling of the Landau bands. 

MSC 2010: 81Q15; 81Q20; 81V70 ; 81S05. 
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1 Introduction 

Schrodinger operators with periodic potentials and magnetic fields have been fascinating 
physicist and mathematicians for the last decades. Due to the competition between the 
crystal length scale and the magnetic length scale, these operators reveals striking features as 
fractal spectrum |GPPO00] . anomalous localization |OGM08] . quantization of the transverse 
conductance |TKNN82| IASS83| IBSE94] . anomalous thermodynamic phase diagrams |Avr04t 
[QXOI] . 

Due to its relevance for the Quantum Hall Effect (QHE) |Mor88t IGra07) . we focus on 
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the two-dimensional Bloch-Landau Hamiltonian 

Hbl:=^ (-^W. -i.^exArj +Vr(r), (1) 

acting in the Hilbert space "Hphy = L'^{M?^d?r), r = (ri,r2) G M^. Here c is the speed of 
hght, h := 271'H is the Planck constant, m is the mass and q the charge (positive if tg = 1 
or negative if Lg = —1) of the charge carrier, B is the strength of the external uniform 
time-independent magnetic field, e± = (0, 0, 1) is a unit vector orthogonal to the sample, 
and Vr is a periodic potential describing the interaction of the carrier with the ionic cores 
of the crystal. For the sake of a simpler notation, in this introduction we assume that the 
periodicity lattice F is simply Z^. 

While extremely interesting, a direct analysis of the fine properties of the operator Hbl 
is a formidable task. Thus the need to study simpler effective models which capture the 
main features of ([I| in suitable physical regimes, as for example in the limit of weak (resp. 
strong) magnetic field. The relevant dimensionless parameter appearing in the problem is 
hs '■= ^o/z^fl oc B~^, where $o = '^'/^ is the magnetic flux quantum, $5 = Q^B is the flux 
of the external magnetic field through the unit cell of the periodicity lattice F (whose area 
is Qr) and Z = l'?l/e is the magnitude of the charge q of the carrier in units of e (the positron 
charge). It is also useful to introduce the reduced constant Hb '■= '^b/2it. 

In the limit of weak magnetic field, Hb — )■ 00, one expects that the relevant features 
are captured by the well-known Peierls' substitution |Pei33l lHar55t lHof76] . thus yielding to 
consider, for each Bloch band S^, = £^=,,(^1,^2) of interest, the following effective model: in 
the Hilbert space L^(T^, d'^k), k being the Bloch momentum and T^ ~ [0, 2tt) x [0, 27r), one 
considers the Hamiltonian operator 

i^Hof ^ =£jk- (^] ^ ex A ^Vk^ ^, V^ G L\T^ dk). (2) 

In physicists' words, the above Hamiltonian corresponds to replace the variables ki and /c2 
in £ri with the symmetric operators {magnetic momenta) 

Since [J^i,J^2\ ^ the latter prescription is formal and ^ must be defined by an appropriate 
variant of the Weyl quantization. The rigorous justification of the Peierls' substitution and 
the definition and the derivation of the Hamiltonian ^ are the content of Section Isl 

The use of the operator ^ traces back to the pioneering works of R. Peierls |Pei33 j and 
P. G. Harper |Har55| . and its spectrum was extensively studied by D. Hofstadter in the 
celebrated paper |Hof 76) . where he specialized to the case £^..{ki, k2) = 2 cos ki + 2 cos A;2- In 
view of that, we call Hofstadter-like Hamiltonian any operator in the form ([2]), while 
the name Hofstadter Hamiltonian is used only for the special case above. However, this 
nomenclature is far to be unique. For instance M. A. Shubin in |Shu94] names discrete 
magnetic Laplacian the same operator (up to a Fourier transform). 
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As for the case of a strong magnetic field, Kb -^ 0, the periodic potential can be con- 
sidered a small perturbation of the Landau Hamiltonian, which provides the leading order 
approximation of Hbl- To the next order of accuracy in Hb-, to each Landau level there 
corresponds an effective Hamiltonian, acting in L^(]R, (ix), given (up to a suitable rescaling 
of the energy scale) by 

i^HarV' = Vrf-«(i,^B)^,xj V^, V'GL'(M,rfx), (4) 

where the r.h.s refers to the ordinary (ig^B)-Weyl quantization of the Z^-periodic function 
Vr : M^ — 7- M. We refer to Section for the definition of the effective Hamiltonian Q. 
We call Harper-like Hamiltonian any operator of the form Q, using the name Harper 
Hamiltonian for the special case Vr(j9, x) = 2cos(27rp) + 2cos(27ra;). 

Remark 1.1 (nomenclature and historical overview). In a remarkable series of papers 
|HS88l IHS90t IHS8 9a] B. Helffer and J. Sjostrand studied the relation between the spectrum 
of the operator dll) and the spectra of a one-parameter family of one-dimensional operators 
on £2(z) defined by 

{hpu)n := Un-i + Un+i + 2 cos(27r6'ra + /3)u„, (5) 

where ^ G M is a fixed number {deformation parameter) and (5 G [0, 27r) is the parameter of 
the family. In the work of the French authors the operator defined by (IS]) is called Harper 
operator (and indeed it was introduced by Harper in |Har55j ). However, in the last three 
decades, the operator ([s]) has been extensively studied by many authors (see |Las05t ILas94) 
for an updated review) with the name of almost- Mathieu operator. To avoid confusion and 
make the nomenclature clear, we chose to adhere to the most recent convention, using the 
name almost- Mathieu operator for ([s]). We thus decided to give credits to Harper's work by 
associating his name to the operator ^. ♦(> 

The regime of strong magnetic field was originally investigated by A. Rauh |Rau74t 
IRau75| . However the correct effective model, the operator (|4]), was derived firstly by M. 
Wilkinson |Wil87] and then, rigorously, by J. Bellissard in an algebraic context |Bel88| and 
by B. Helffer and J. Sjostrand in |HS89bj . inspired by the latter paper. In both these papers 
the case of weak magnetic field is also considered. In particular, in |HS89bj it is proven that 
-f^Har (resp. -^Hof) has, locally on the energy axis, the same spectrum and the same density 
of states of Hbl, as ^^ — ;■ (resp. fiB — !■ oo). 

Beyond the spectrum and the density of states, there are other mathematical prop- 
erties of Hbl which reveal interesting physics, as for example the orbital magnetization 



[GA031 ITCVR051 ICTVR06] or the transverse (Hall) conductance "^^^ These properties are 
not invariant under a loose equivalence relation as isospectrality, then it is important to 
show that -f^Har (resp. i^Hof) is approximately unitarily equivalent to Hbl in the appropriate 



^^•'Some authors refer to the "Hall conductivity" rather than to the "Hall conductance". Indeed, in two space 
dimensions the conductivity (microscopic quantity) coincides exactly with the conductance (macroscopic 
quantity), so the two concepts are synonymous. In this sense the quantization of the Hall conductivity is a 
macroscopic quantum phenomenon. 
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limit. This is the main goal of this paper. The problem is not purely academic, since there 
exist examples of two isospectral Schrodinger operators which are however not unitarily 
equivalent and which exhibit different values of the transverse conductance. In particular, 
in [DFPlOb] we prove that the Hofstadter operator and the Harper operator are isospectral 
but not unitarily equivalent. One concludes that, in the study of the conductance, it is not 
enough to prove that the effective model is isospectral to the original Hamiltonian. 

We thus introduce the stronger notion of unitarily effective model, referring to the 
concept of almost-invariant subspace introduced by G. Nenciu |Nen02] and to the related 
notion of effective Hamiltonian |PST03b] IPST03a[ ITeu03) . which we shortly review. Let us 
focus on the regime of weak (resp. strong) magnetic field and define e = e{B) = ^/tiB (resp. 
e{B) = Kb) so that e — )■ in the relevant limit. Let H^ be an orthogonal projector in "Hphy 
such that, for any N ENyN < Nq there exist a constant Cn such that 

||[HBL;n,]||<C;v£'^ (6) 

for e sufficiently small. Then RanH^ is called an almost-invariant subspace |Nen02t 
ITeuOS] at accuracy A^o, since it follows by a Duhammel's argument that 



H,) e-^^"«^ HJI <CNe' 



for every s G M, A^ < A^o- Granted the existence of such a subspace, we call (unitarily) 
effective Hamiltonian a self-adjoint operator ifgfr acting on a Hilbert space "Href? such 
that there exists a unitary f/^ : RanH^ — )■ "Href such that for any A^ G N, A^ < A"o, one has 

II (H, Hbl - U-'H,sU,) H,|| < C^ e^. (7) 

The estimates (|6| and ([T]) imply that 

II ^g-»sHBL _ U-^e-''^^«Ue) Hell < Cl^ e^\s\. (8) 

When the macroscopic time-scale t = es is physically relevant, the estimate above is simply 
rescaled. The triple (T^ref, U^, Hes) is, by definition, a unitarily effective model for Hbl- To 
our purposes, it is important to notice that the asymptotic unitary equivalence in ([T]) assures 
that the topological invariants related with the spectral projections of H^ HgL n^ (ii'-theory, 
Chern numbers, . . . ) are equal to those of H^s, for e sufficiently small [DP09] . 

In this paper we prove that in the limit Kb ^ oo the Hofstadter-like Hamiltonian (pi) 
provides a unitarily effective model for Hbl with accuracy A"o = 1, and we exhibit an iterative 



algorithm to construct an effective model at any order of accuracy A'q G N (Theorem 3.12). 
As for the limit Hb -^ C)0, up to a rescaling of the energy, the non-trivial leading order 
(accuracy A'o = 1) for the effective Hamiltonian is given by the Haper-like Hamiltonian Q. 
We also exhibit the effective Hamiltonian with accuracy A^o = 2, i.e. up to errors of order 



0{e^) (Theorem 4.7). Moreover, due to the robustness of the adiabatic techniques, we can 
generalize the simple model described by ([I| to include other potentials, see (|9]), including 
in particular a periodic vector potential Ap. This terms produces interesting consequences 



especially in the Harper regime (see Section 4.8) and it could play a relevant role in the 
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theory of orbital magnetization. This kind of generahzation is new with respect to both 
[BilSS] and |HS89b| . 

Our proof is based on the observation that both the Hofstadter and the Harper regime are 
space-adiabatic hmits, and can be treated in the framework of space- adiabatic perturbation 
theory (SAPT) |PST03b( IPST03a) . see also |Teu03] . As for the Hofstadter regime, the 
proof follows ideas similar to the ones in |PST03a] . Our generalization allows however to 
consider a constant magnetic field (while in |PST03a] the vector potential is assumed in 
C^(]R'^)) and to include a periodic vector potential. Moreover the proof extends the one 



in |PST03a] . in view of the use of the special symbol classes defined in Section 3.4 On 



the contrary, from the discussion of the Harper regime Hb ^ some new mathematical 



problems emerge. Then, although the "philosophy" of the proof of Theorem 4.7 is of SAPT- 
type, the technical part is new as it will be explained in Section |4j Notice that the regime 
of weak magnetic field can also be conveniently approached by using the magnetic Weyl 
quantization |MP04l IMPR051 IIMP07lllMP09] . a viewpoint which is investigated in |DL10] . 

For the sake of completness, we summarize some salient aspects of the SAPT. We refer 
to specific references (e.g. |Teu03) ) for a complete exposition. Let H be the Hamiltonian 
of a generic physical system which acts on the total (or physical) Hilbert space "H. For the 
SAPT to be applicable, three important ingredients needs: (i) a distinction between fast 
and slow degrees of freedom which is mathematically expressed by a unitary decomposition 
of the physical space 1-i into a product space T-Ls ® ^f (or, more generally, a direct integral), 
the first factor being the space of the slow degrees of freedom and the second the space of 
the fast degrees of freedom; Tis = L'^{A4) for suitable measure space Ai is also required; (ii) 
a dimensionless adiabatic parameter e -C 1 that quantifies the separation of scales between 
the fast and slow degrees of freedom and which measures how far are the slow degrees of 
freedom to be "classical" in terms of some process of quantization; (iii) a relevant part of 
the spectrum for the fast dynamics which remains separated from the rest of the spectrum 
under the perturbation caused by the slow degrees of freedom. 

A numerical simulation of the spectrum of the Hofstadter operator, as a function of the 
parameter e = Y^b, leads to a fascinating picture known as Hofstadter butterfly |Hof 76] . 
Since the spectrum has zero measure as a subset of the square, the physically relevant object 
is its complement, the resolvent set. It has been pointed out by D. Osadchy and J. Avron 
[OAOlj that the open connected regions of the resolvent set (islands) can be associated to 
different thermodynamic phases (at zero temperature) of a gas of non interacting fermions in 
a periodic potential, with e oc B and the chemical potential as thermodynamic coordinates. 
The different phases are labeled by an integer (topological quantum number), interpreted 
as the value of the transverse conductance of the system in units of <^^/h in the limit of 
weak magnetic field. The latter integers are conveniently visualized by different colors, thus 
leading to the colored Hofstadter butterfly |OA01l IAvr04| . With this language in mind, 
the main result of this paper can be reformulated by saying that the Hofstadter-like and 
Harper-like Hamiltonians are "colour-preserving effective models" for the original Bloch- 
Landau Hamiltonian. Thus they describe, though in a distorted and approximated way, 
some aspects of the thermodynamics of the original system. 

Acknowledgments. It is a pleasure to thank Y. Avron, J. Bellissard, G. Dell'Antonio, 
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2 Description of the model 

A generalized Bloch-Landau Hamiltonian 

The Hamiltonian ([I| describes the dynamics of particle with mass m and charge q which in- 
teracts with the ionic structure of a two dimensional crystal and with an external orthogonal 
uniform magnetic field. A more general model is provided by the operator 



2m 



1 r . — ?« / X Q ^^ 



-ihVr - -Ar (r) - -A (r) 

c c 



+ Vr (r) + g $ (r) (9) 



still called Bloch-Landau Hamiltonian and, with an abuse of notation, still denoted with 
the same symbol used in ([I|. The vector-valued function A := (Ai, A2) is a vector potential 
corresponding to an (orthogonal) external magnetic field B = V^ A A = (9iA2 — 92A1) e_L, 
$ is a scalar potential corresponding to a (parallel) external electric field E = — V^^ and 
Ar and Vr are internal periodic potentials which describe the electromagnetic interaction 
with the ionic cores of the crystal lattice. The external vector potential is assumed to have 
the following structure 

A(r)=Ao(r) + AB(r), (10) 

where Aq is a bounded function and A^ describes a uniform orthogonal magnetic field of 
strenght B, i.e. in the symmetric gauge 

B / B B \ 

Aij(r) = -exAr= (--r2,-rij , V, A Ab = 5 e^, V, ■ A^ = 0. (11) 

The evolution of the system is prescribed by the Schrodinger equation 

i^— ^(r,s) = HBL V'(r,s), (12) 

where s corresponds to the microscopical time-scale. 

Mathematical description of the crystal structure 

The periodicity of the crystal is described by a two dimensional lattice F C M^ (namely a 
discrete subgroup of maximal dimension of the Abelian group (]R^,+)), thus F c^ 1? . Let 
{a, 6} C M^ be two generators of F, i.e. 

F = {7 G M^ : 7 = n\a + ^26, "n-i, n^ G Z}. 

The fundamental or Voronoi cell of F is M-p := {r & M.'^ \ r = li a + I2 b, li,l2 G [0, 1]} 
and its area is given by fir = |a A 6|. We fix the orientation of the lattice in such a way that 
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VLy = {aib2 — a2bi) > 0. We say that a function /r : M^ — )■ C is T-periodic if fr{f' + 'y) = fr{f) 
for all 7 G r and all r G M^. The electrostatic and magnetostatic crystal potentials Vr and 
Ay are assumed to be F-periodic according to the previous definition. 

An important notion is that of dual lattice T* which is the set of the vectors 7* G M^ such 
that 7* -7 G 27rZ for all 7 G T. Let {a*, 6*} C M^ be defined by the relations a*-a = b*-b=l 
and a* ■ b = b* ■ a = 0; these vectors are the generators of the lattice F*, i.e. 

r* = {7* eR"^ : 7* = mi 27ra* + ma 27Tb*, mi, ma G Z}. 

The Brillouin zone Mr* := {k &M? \k = ki a* + ^2 6*, /ci, ^2 G [0, 27r]} is the fundamental 
cell of the dual lattice T* . The explicit expressions for the dual generators {a*, b*} in terms 
of the basis {a, 6} is 

|a A o| Uy \a/\b\ \Iy 



It follows from (13) that the surface of the Brillouin zone is Vt^* = (27r) |a* /\b*\ = (2^) /n 



Given a F-periodic function /r, we denote its Fourier decomposition as 

/r(r) = Y. /(^*) ^*'*" = E ^-1.™^ ^''"^'"' "*^'"' '*^"'- (14) 

7*£r* mi,m2SZ 

A 1? -periodic function / : M^ ^^ C is a function periodic with respect to an orthonormal 
lattice, namely such that /(xi + 1, X2) = /(xi, 0:2 + 1) = /(xi, 0:2) for all Xi, X2 G M. If one 
replaces the two real variables by Xi := a*-r and 0:2 := b*-r one has that /r(^) := f{o*-^i b*-r) 
is F-periodic in r. Every F-periodic function can be obtained in this way. 

Assumptions on the regularity of the potentials and self-adjointness 

Let us denote by C^(]R^, M) the space of real- valued n-times differentiable functions (smooth 
functions if n = 00) with continuous and bounded derivatives up to order n. Concerning 
the internal potentials Ar and Vr we need to assume that: 

Assumption (As) [internal potentials, strong form]. The T-periodic potential \/r and 
the two components of the T-periodic vector potential Ap are functions of class C^(]R^,]R). 

Sometime will be enough to consider a weaker version of this assumption, namely: 

Assumption (Aw) [internal potentials, weak form]. The two components of the 
T-periodic vector potential Ar are in Cl(R^,M.). The T-periodic potential Vr verifies the 
condition J^ |Vr(r)p d'^r < +00. 

Assumption (A^) implies that Vr is uniformly locally L^ and this implies also that Vp 
is infinitesimally bounded with respect to — A^ (see Theorem XIII. 96 in |RS78| ). Concern- 
ing the external potentials A and $, we need to assume that: 
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Assumption (B) [external potentials]. The scalar potential $ is of class C^ 



The vector potential A consists of a linear term Ab of the form (11) plus a bounded term 
Ao which IS of class ^^(M^M). 

When the external potentials A and $ vanish, the Bloch-Landau Hamiltonian ^ reduces 
to the periodic Hamiltonian (or Bloch Hamiltonian) 



H 



per • 



2m . 



-iKVr 



^Ar(r) 



+ Vr(r) 



(15) 



The domains of self-adjointness of Hbl and Hper are described in the following proposition. 
Its proof, together with some basic notion about the Sobolev space [K^(]R^) and the magnetic- 
Sobolev space CK^(]I 



is postponed to Section A.l 



Proposition 2.1. Let Assumptions (A^) and (B) hold true. Then both Hbl and Hper 
are essentially self-adjoint operators on L^(]R^,(i^r) with common domain of essential self- 
adjointness the space of smooth functions with compact support C^(]R^,C). Moreover the 
domain of self-adjointness of Hper is 'K'^(M?) while the domain of self-adjointness of Hbl 

zsM^(M2). 



3 Space-adiabatic theory for the Hofstadter regime 
3.1 Adiabatic parameter for weak magnetic fields 



The SAPT for a Bloch electron developed in |PST03a] is based on the existence of a sep- 
aration between the microscopic space scale fixed by the lattice spacing i := a/Ot, and 
a macroscopic space scale fixed by the scale of variation of the "slowly varying" external 
potentials. The existence of such a separation of scales is expressed by introducing a di- 
mensionless parameter £ ^ 1 {adiabatic parameter) to control the scale of variation of the 
vector potential and the scalar potential $ appearing in (pi), namely by setting A = A (err) 
and $ = $(£r). In particular the external magnetic and electric fields are weak compared 
to the fields generated by the ionic cores. 

It is useful to rewrite the (e-dependent) Hamiltonian (|9| in a dimensionless form. The 
microscopic unit of length being i, we introduce the dimensionless position vector x := ^/e 
and the dimensionless gradient V^ = iVr- Moreover, since the vector potential has the 
dimension of a length times a magnetic field, then A{ex) := "^/m A{eix) is a dimensionless 
function, with B a dimensional constant which fixes the order of magnitude of the magnetic 
field due to the external vector potential A. Similarly for Ap (with e = 1). Factoring out 
the dimensional constants one finds 



Hbl '■— t^Hbl 

•^0 



1 
2 



-1 2 



-iVx 



qnrBi 




Ar (x) - Lq ^ A {ex) 

en ,e 

--n- 



+ rr(x) + 0(ex), (16) 
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where £q := ^^/mnr is the natural unit of the energy fixed by the problem, Vy{x) := ^/sq Vr(£x) 
and ^(ex) := i/so ^{elx) are both dimensionless quantities. The constant Hy will play 
no particular role in the rest of this paper, so it is reabsorbed into the definition of the 
dimensionless vector potential Ap, i.e. formally Hy = 1. 

Comparing the dimensional Hamiltonian (16) with the original Hamiltonian ^ , or 



observing that the strenght of the magnetic field goes to zero (at least linearly) with e, it is 
physically reasonable to estimate eHb oc 1. This is rigorously true in the case in which the 
external electromagnetic field is uniform. 

The external force due to A and (p are of order of e and therefore have to act over a 
time of order e^^ to produce a finite change, which defines the macroscopic time-scale. The 
macroscopic dimensionless (slow) time-scale is fixed hj t := e^s where s is the dimensional 



microscopic (fast) time-scale. With this change of scale the Schrodinger equation (12) reads 



iSj^iP = Hbl ^ (17) 



with i^BL given by equation (16). 



Remark 3.1. Observe that from the definition of the dimensionless periodic potential Ar 
and Vr it follows that they are periodic with respect to the transform a; i— )■ a; + i/i. This 
means that A^ and Vp are periodic with respect to a "normalized" lattice whose fundamental 
cell has surface 1. ♦() 

3.2 Separation of scales: the Bloch-Floquet transform 

To make explicit the presence of the linear term of the external vector potential, we can 



rewrite the (16) as follows 



i^BL=2 



-zV — Ar {x) — Aq {ex) — tg -e± A ex 



+ Vr (x) + {ex) , (18) 



where the adiabatic parameter e expresses the separation between the macroscopic length- 
scale, defined by the external potentials, and the microscopic length-scale, defined by the 
internal F-periodic potentials. The separation between slow and fast degrees of freedom can 
be expressed decomposing the physical Hilbert space "Hphy = L^{M.'^,d'^x) into a product 
of two Hilbert spaces or, more generally, into a direct integral. To this end, we use the 
Bloch-Floquet transform |Kuc93) . As in |PST03a] we define the (modified) Bloch-Floquet 
transform Z oi a function ip G 5(]R^) to be 

(ZV^)(A;,0) :=^e^*(^+^)-V(^ + 7), {k,e) eR"^ x R"^. (19) 

7er 

Directly from the definition one can check the following periodicity properties: 

{Zij){k,9 + j) = {Zij){k,9) V7Gr (20) 

{Zilj){k + ^*,e)=e-''-^*{Zilj){k,e) V7*eF*. (21) 

10 
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Equation (20) shows that for any fixed k E 



M2, {Z^){k, 
with V := K 



■) is a F-periodic function and can 
^/r a two-dimensional slant torus 



be seen as an element of "Hf 
( Voronoi torus). The torus V coincides with the the fundamental cell Mp endowed with the 
identification of the opposite edges and (P6 denotes the (normalized) measure induced on V 
by the identification with M-p. The Hilbert space Tif is the space of fast degrees of freedom, 



corresponding to the microscopic scale. Equation (21) involves a unitary representation 
T : r* — > ^{l-ii) of the group of the (dual) lattice translations F* on the Hilbert space 
"Hf. For every 7* G F* the unitary operator r(7*) is the multiplication with e^^"^ 
convenient to introduce the Hilbert space 



It will be 



rL-T 



{^eLl^{R\d%n,) : ^{k-^\-) = r{Y)Hk,-)} 



(22) 



d'^k 



equipped with the inner product {ip;ip)-}i^ := Jj^^ ^{ilj(k);if{k))fi^ d^k where d^k :— j^^ 

is the normalized measure. There is a natural isomorphism from H^ to L^ (Mr*,(i^A;, Hf) 
given by restriction from M? to Mp* , and with inverse given by r-covariant continuation, as 



suggested by (21). The Bloch-Floquet transform ([I9j) extends to a unitary map 



I-phy 



(23) 



The Hilbert space L"^ {Mr* , d"^ k) can be seen as the space of slow degrees of freedom and in 
this sense the transform Z produces a decomposition of the physical Hilbert space according 
to the existence of fast and slow degrees of freedom. 

We need to discuss how differential and multiplication operators behave under Z. Let 
Q = {Q11Q2) be the multiplication by a: = (xi,X2) defined on its maximal domain and 
—iVx with domain the Sobolev space J{^{I 



P={Pl,P2 



Z P z 



-1 



k(g)l'Uf + lL2(A^r.) ® -'i'^e, 



, then from (19) it follows: 
Z Q Z-^ 



^vi 



(24) 



where —iVg acts on the domain !K^(V) while the domain of the differential operator iV^ 
is the space Tir H JCj^^ (M^j'Hf), namely it consists of vector- valued distributions which are 



in "K^ (Mrtj'Hf) and satisfy the ^-dependent boundary condition associated with (21). The 



central feature of the Bloch-Floquet transform is, however, that multiplication operators 
corresponding to F-periodic functions like A^ or Vp are mapped into multiplication operators 
corresponding to the same function, i.e. 



Z Ar{x) Z-^ = lL2(^Mr-.) ® M0) 



Z Vr(x) Z- 



1 



L2 



(M,,)®Vr{9). (25) 



Let H^ := Z H^i^ Z ^ he the Bloch-Floquet transform of the Bloch-Landau Hamiltonian 



(18). According to the relations (24) and (25) one obtains from (18) that 

1 



H' 



-iVe + k~Ar{e)- A^ {leVD 



:e± A {teVl] 



+ Vr{e) + <p{teVl) (26) 



with domain of self-adjointness Z'K^^i 
magnetic-Sobolev space. 



^) C Tir, i.e. the image under Z of the second 
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3.3 The periodic Hamiltonian and the gap condition 



When e = the Bloch-Landau Hamiltonian (18) reduces to the periodic Hamiltonian 

^per = \ [-^V,. - Ar {x)f + Vr (x) . (27) 



According to (26) the Bloch-Floquet transform maps -ffper into a fibered operator. In other 
words, denoting if^j. := Z Hper Z~^, one has H^^^ = J^ ^ H^erik) d?k where, for each 
k e My* 

Hperik) = ^ [-iVe + k- Ar{e)f + Vr{e). (28) 

The operator H^evik) acts on 7{f = L'^{Y,(P6) with self-adjointness domain V := J{^(V) 
(the second Sobolev space) independent of A; G Mr* . Moreover it is easy to check that the 
Bloch-Floquet transform induces the following property of periodicity, called T-equivariance: 

Hperilk] - 7*) = r(7*) Hperilk]) t{^T' ^ T* V7* G r . (29) 

where the notation k := [k] — 7* denotes the a.e.-unique decomposition of A; G M^ as a sum 

of [k] G Mr* and 7* G T*. 

Remark 3.2 ( Analiticity) . For any k G M^, let I(k) be the unitary operator acting on 
"Hf as the multiplication by e~*^'^. Obviously I{k) = I{[k] — 7*) = /([A;])r(7*)~^ A simple 
computation shows that 

Hperik) = I{k) HpeM I{k)-' (30) 

where the equality holds on the fixed domain of self-adjointness V = J{^(V). The r- 



equivariance property (29) follows immediately from (30). Moreover from (30) is evident 
that iJper(fc) defines an analytic family (of type A) in the sense of Kato (see |RS78] Chapter 
XII). Finally a short computation shows 



(4,/^per)(A:) = -tl{k) [df /7per(0)] I{k)-^ = I{k) i-zVe - Ar (6))^ I{k)-' 
and {dlMper){k) = Iv, {dl^^^Hp^r){k) = on the domain V. ♦<> 

The spectrum of -f^per, which coincides with the spectrum of if£,j., is given by the union 
of all the spectra of Hper{k). The following classical results hold true: 

Proposition 3.3. Let Vp and Ar satisfy Assumption (A„), then: 



(i) for all A; G M^ the operator Hper{k) defined by (27) is self-adjoint with domain V = 
[K^(V) and is bounded below; 

(ii) Hper{k) has compact resolvent and its spectrum is purely discrete with eigenvalues 
^nik) — !■ +00 as n ^ +00; 

(iii) let the eigenvalues be arranged in increasing order and repeated according to their 
multiplicity for any k G M^*, i.e. £i{k) ^ £2{k) ^ £z{k) ^ ... then Sn{k) is a 
continuous T* -periodic function of k. 
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The above result differs from the standard theory of periodic Schrodinger operators just 
for the presence of a periodic vector potential Ay- Since we were no able to find a suitable 
reference in the literature, we sketch its proof in Appendix A.l 



We call £n{-) the n^^ Block hand or energy hand. The corresponding normalized eigen- 
states {v5n(^)}neN C T> are called Block functions and form, for any k G Mr*, an orthonormal 
basis of "Hf. Notice that, with this choice of the labelling, £n{') and <^„(-) are continuous in 
/c, but generally they are not smooth functions if eigenvalue crossings are present. 

We say that a family of Bloch bands {^n(')}n6X) with X := [/+, /_] fl N, is isolated if 
Jnf dist [ \J{£n{k)}, \J{m)} I = C, > 0. (31) 

The existence of an isolated part of the spectrum is a necessary ingredient for an adiabatic 
theory. We introduce the following: 

Assumption (C) [constant gap condition]. Tke spectrum of Hp^^ admits a family 



of Block hands {£n{')}nex wkick is isolated in tke sense of (31). 



Let Pi{k) be the spectral projector of ifper(/c) corresponding to the family of eigenval- 
ues {Sn{k)}n(^x, then P-f '■= j\j ^ Px{k) d'^k is the projector on the isolated family of Bloch 
bands labeled by X. In terms of Bloch functions (using the Dirac notation), one has that 
Px{-) = Ylin(^x\Vn{.-)){Vn{.-)\- Howcvcr, in general, (^n(-) are not smooth functions of k at 
eigenvalue crossing, while Px(-) is a smooth function of k because of the gap condition. 
Moreover, from the periodicity of ifper(-), one argues Px{[k] — 7*) = t(7*) Px{[k]) r(7*)~^. 
In general the smoothness of Px(-) is not enough to assure the existence of family of or- 
thonormal basis for the subspaces RanPi(-) which varies smoothly (or only continuously) 
with respect to /c G Mp*. Then we need the following assumption. 

Assumption (D) [smooth frame]. Let {£n{-)}nex he a family of Block hands (\X\ = 
m > 1). We assume tkat tkere exists an ortkonormal hasis {ipn{-)}lLi o/ RanPj(-) wkose 
elements are smootk and (left) r-covariant witk respect to k, i.e. iIjj{- — 7*) = t(7*)V^j(-) for 
all j = 1,. . . ,m and 7* G F* . 

Note that it is not required that ipj{k) is an eigenfunction of Hperik). However, in the 
special but important case in which the family of bands consist of a single isolated m-fold 
degenerate eigenvalue, i.e. Sn{k) = £*{k) for every n = 1, . . . ,m, then the Assumption (D) 
is equivalent to the existence of an orthonormal basis consisting of smooth and r-covariant 
Bloch functions. 

Remark 3.4 (Time-reversal symmetry breaking). As far as low dimensional models are 
concerned ((i < 3), Theorem 1 in |Pan07| assures that Assumption (D) is true whenever the 
Hamiltonian Hp^r is invariant with respect to the time-reversal symmetry., which is imple- 
mented in the Schrodinger representation by the complex conjugation operator. However, 
the term Ar 7^ in Hp^r generically breaks the time reversal symmetry. Therefore, to con- 
sider also the effects due to a periodic vector potential, we need to assume the existence of 
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a smooth family of frames. Anyway is opinion of the authors that the result in |Pan07| can 
be extended to the case of a periodic vector potential, at least assuming that A^ is small in 
a suitable sense. ♦(> 

Let ko be a fixed point in Mr* and define the projection tt^. := Px{ko). If the Assumption 
(C) holds true then dim tTj. = dim Px{k) for all k G M^. Let {Xn}T=i be an orthonormal 
basis for Ran VTr and define a unitary map 

Uo{k) ■=UQ{k) + UQ{k), with Uo{k) := ^ \Xj){^j{k)\, (32) 

i<j<e 

which maps Ran Px{k) in Ran VTr. The definition of this unitary is not unique because 
the freedom in the choice of the frame and of the orthogonal complement UQ{k). From 
the definition and the r-covariance of ipj{-) one has that Mo(fc) Px{k) UQ{k)~^ = tt,. and 
Uo{[k] — 7*) = Uo{[k]) r(7*)^^ {right r-covariance). 

3.4 r-equivariant and special r-equivariant symbol classes 



Proposition 3.3 shows that for all k G M^, the operator Hpej.{k) defines an unbounded self- 
adjoint operator on the Hilbert space "Hf with dense domain V := !K^(V). However the 
domain P can be considered itself as a Hilbert space with respect to the Sobolev norm 
II ■ \\v '■= ||(lwf — Ag) ■ ll^j and so H^er{k) can be seen as a bounded linear operator from 
T) to Hf, i.e. as an element of the Banach space ^{Vj'Hf). The map M^ 9 /c h- )■ H^erik) G 
SS{T>^ T-Lf) is a special example of a operator-valued symbol. For a summary about the theory 
of the Weyl quantization of vector-valued symbols, we refer to Appendices A and B in 



|Teu03 ] . In what follows we will need the following definition. 

Definition 3.5 (Hormander symbol classes). A symbol is any map F from the (cotangent) 
space ]R2 X M^ ^ the Banach space ^{V,ni), i.e. M^ x ]R2 3 (^^^) ^ F{k,r]) G ^{V^Ui). 
A function u; : M^ x M^ — >• [0, +00) is said to he an order function if there exists constants 
Co > and Nq> such that 



w 



{k, r])^Co{l + \k-k'\^ + \7]- r^f ) ^ w{k', 7]') (33) 



for every {k,r]), {k',r]') G M^ x M^ j[ symbol F G C°°(]R2 x R'^,^{V,H{)) is an element of 

the (Hormander) symbol class S'^{^{V,'Hf)) with order function w, if for every a,/3 G N^ 

there exists a constant Ca,j3 > such that \\{d'j^d!^F){k,rj)\\^(x>,Hi) ^ Ca,i3w{k,ri) for every 
(A;,r7) gM^ xM2_ 

According to the previous definition, the vector- valued map ifper(-) defines a Hormander 
symbol constant in the r^- variables and with order function v{k, rj) := 1 + 1 A;p (see the proof of 



Proposition 3.10 below). However, as showed by equation (29), the symbol i7per(-) satisfies 



an extra condition of periodicity. 

Definition 3.6 (r-equivariant symbols). Let T* be a two dimensional lattice (the dual 
lattice defined in Section [^ for our aims) and r : T* ^ '^ (Hi) the unitary representation 
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defined in Section 3.^ Denote by r := r|p the bounded- operator ^ ' representation ofT* in 



T>. A symbol F G S'^{^{V,'H{)) is said to be r-equivariant if 

F{k - 7*, r/) = r(7*) F{k, 7]) r(7*)"^ V 7* G P, A; G Ml 

The space of r-equivariant symbols is denoted as S^{^{V, !-[{)). 

For the purposes of this work, it is convenient to focus on special classes of symbols. By 
considering the kinetic momentum function M^ x M^ 9 (/c,?]) 1 — > K,{k,ri) := k — A{r]) G M^, 
with A fulfilling Assumption (B), one defines the minimal coupling map by 

{k, v) ^ Jn{k, T]) := {K{k, r/), r/) G M^ X M^. (34) 

Definition 3.7 (Special r-equivariant symbols). Let w be an order function, in the sense 



of (^. We define 

Sl^{^{V,n,)):={F = Foj, : F e S^^i^iV.U,))}. 

We refer to S^.^{^{T>,'H^)) as the class of special r-equivariant symbols. The following 
result shows that special symbols can be considered as genuine r-equivariant symbols with 
respect to a modified order function. The key ingredient is the linear growth of the kinetic 
momentum. 

Lemma 3.8. With the above notations S^.^{^ {V , rCf)) C Sf{^{V,nf)) where w' := woj^. 

Proof. If F G ^^(.^(PjT/f)) then also Foj^ is r-equivariant, indeed K,{k—'~f*, t]) = K,{k, t]) — 
7* and (F o jK){k — 'j*,r]) = r{Y) (-^ ° jK){k,r]) r{'-y*)~'^. Since j^ is a smooth function, 
then also the composition F o j^ is a smooth function. Observing that {F o jK){k,f]) = 
F{k — A{'r]),ri) it follows that 

{d,^{Foj^))(k,r]) = i{d,^F)oj^){k,r^) 

2 
{^,^{FoJ,)){k,v) = {{^,^F)oJ^){k,r^) + J2i^,^^)ik,v){i^k.F)oJ,){k,v) 

i=\ 

where d^.Ki are bounded functions because Assumption (B). From the first equation it 
follows that 

||5fc,(i^ojK)(A;,^)||,«(©,Wf) ^ ^,0 iwoj^){k,r]) j = 1,2 

for suitable positive constants Cjo- Similarly the second equation implies 

\\drj,{F o j^){k,ri)\\.j(^v,nf) < [Coj + K{Ci^o + C2,o)]{w o j^){k,ri). 

where fi' > is a bound for the functions dr^Ki- By an inductive argument on the number 
of the derivatives one can proof that the derivatives of F o j^ are bounded by the function 
w' := w o ji^. To complete the proof we need to show that w' is an order function according 



to Definition 3.5 This follows by a simple computation using the fact that k has a linear 



growth in k and r]. 



•^^'Clearly r(7*) acts as an invertible bounded operator on the space V. but it is no longer unitary with 
respect to the Sobolev-norm defined on T). 
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In view of Lemma 3.8 all the results of Appendix B of |Teu03| hold true for symbols in 
S^.^{^{T>,'Hf)) and in particular the quantization of a symbol in 5'^^(=^(X','Hf)) preserves 
the r-equivariance and that the pointwise product or the Moyal product of two symbols of 
order Wi and W2 produce a symbol of order W1W2 (see |Teu03) . Propositions B.3 and B.4). 

Remark 3.9 (Notation). In what follows we will use the short notation F{n\ri) := {F o 
jK){k,ri) to denote the special symbol F o j^ g S'^.^{^(V,'Hf)) related to the r-equivariant 
symbols F G S'^(=^(P,'Hf)). We emphasize on the use of the semicolon ";" instead the 
comma "," and of the symbol of the kinetic momentum k instead the quasi-momentum k. 

♦0 



3.5 Semiclassics: quantization of equivariant symbols 



As explained in Section 3^, the Bloch-Floquet transform Z provides the separation between 
the fast degrees of freedom, associated to the Hilbert space "Hf = L'^iY ,d?9), and the slow 
degrees of freedom, associated to the Hilbert L"^ {My* , (P k) . A fruitful point of view is to 
consider the slow degrees of freedom "classical" with respect to the "quantum" fast degrees of 
freedom. Mathematically, this is achived by recognizing that the Hamiltonian H^ defined in 



( 26 ) is the Weyl quantization of an operator- valued "semiclassical" symbol over the classical 
phase space M^ x M^. As explained rigorously in the Appendices A and B of |Teu03| . the 
quantization procedure maps an operator-valued symbol F : M^ x M^ — >■ ^(P, "Hf) into a 
linear operator Op^(F) : S(M?,V) -^ S(M?,ni), where S{R'^,n) denotes the space of Ti- 
valued Schwartz functions. The quantization procedure concerns only the slow degrees of 
freedom and at a formal level can be identified with the prescription 

k I — )■ Op^{k) := multiplication by A; (g) Ix?; r] 1 — )■ Op^{ri) := ieVk ® tv- (35) 

Let us consider the operator-valued symbol Hq : M."^ x M.^ ^ ^{V, Tii) defined by 



Hoik,v):=l 



1 ^' 

-iVe + k-Ar{e)- Aq^t]) - Lg-e± A r/ 



+ Vri9) + (P{r]). (36) 



The symbol Hq does not depend on e and in view of Proposition 3.3 it defines an unbounded 



operator on T-Lf with domain of self-adjointness V = "K (V) for all choice of (k, rj) G 



According to the notation of Section 3.4 and comparing ^36^ with ^28^ we can write 



H,{k, ri) = Hp,, {Kik, r])) + ^(ri) = {H^ o j,)(A;, r/). (37) 



where H^{k,ri) := H^erik) + (pirf). As suggested by equation (29), iJ<^ is a r-equivariant 
symbol. Thus the symbol Hq is r-equivariant with respect to the kinetic momentum k. The 
following result establishes the exact symbol class for Hq. 

Proposition 3.10. //Assumption (A^) and (B) hold true thenH^ G S^.^{^(V;'Hi)) with 
order function v{k,ri) := 1 + |A;p. 



Proof. Using the result of Lemma 3.8, we only need to show that H^ G S'"{J^{V, Hi)). The 



later claim is easy to verify, indeed the derivative in tj are bounded functions, the second 
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derivative in A; is a constant and the derivatives of higher order in k are zero. Then we need 
only to check the growth of the first derivative in k. A simple computation shows that 



'fc,^<^)(^,^)l|.^(2?,«f) 



^H^er){k) {ln,-^t 



\-ii 



\.n'Hi) 



and since d^ Hp^r is r-equivariant (see Remark 3.2), then 



imHt)- 



Observing that r(7*) is the multiplication by e*^'''* in Tif and by a simple computation that 
{dk^Hper){[kMY)-' = r(7*)-i[-27; + {dk^Hp,,){[k])] one has 



where Ci = 2||(l^j-Ae) ^\\fij{nf), C2 ■= max^g 



Mr 



ij -n per 



,«(D,-Hf) 



and |7*| < I7* 



1^ — [^]| ^ 1^1 + C's with C3 := maxfcgAf^, \k\. The claim follows observing that 1 + |/i;| ^ 
2(l + |A;p). ■ 



a(Ho) 




Figure 1: Structure of the spectrum of Ho{k, 77). The picture shows schematically a "relevant part of the 
spectrum", consisting of two energy bands {E^,, i?*+i}, with E^,^j{k, -q) = £^,+j{K{k, 77)) + (l>{rj). Notice that 
we assume only a local gap condition, as stated in pS] ), while in the picture a stronger condition is satisfied: 
a gap exists when projecting the relevant bands on the vertical axis. 



Equation (37) provides information about the dependence on k and i] of the spectrum 



of Hq. The n eigenvalue En{k, rj) of the operator Ho(k, rf) is related to the n eigenvalue 
8n{k) of the periodic Hamiltonian Hper{k) by the relation En{k,r]) = £n{n{k,ri)) + (pirj). 



The function E^ '■ 



-)■ 



is still r*-periodic in k but only oscillating with bounded 
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variation in rj. Assumption (C) for the family of Bloch bands {£n{-)}nex immediately implies 
that 

inf dist\[j{E^{k,r^)},[j{E,{k,r^)}\ = C, > 0. (38) 

{fc,r,)eAfp.xM2 \nYx m ) 

This is the relevant part of the spectrum of Hq which we are interested in. 

According to the general theory (see Appendices A and B of |Teu03] ). one has that 

2 

+ VT{e)+(t){ieVk) (39) 



Op.(i^o) = \ 



-tVe + k- Ar{e) - Aq (ieVk) - 6g-e± A (teVk) 



defines a linear operator from cS(M^, D) in 5(M^, T-Lf) and by duality it extends to a continuous 
mapping Op^^^Ho) : 5'(]R^,'D) — )■ S'(M.'^,7if) (with an abuse of notation we use the same 
symbol for the extended operator). The r-equivariance assures that Op^{Ho)(f{[k] — 7*) = 
r(7*)0pg(ifo)v^([^]) (see [Teu03j Proposition B.3). Since Op^{Hq) preserves r-equivariance 
it can then be restricted to an operator on the domain 2^CK^(]R^) C iS'(]R^,X') which is the 



domain of self-adjointness of H^, according to (26). To conclude that Op^{Hq), restricted 
to 2CK^(M^), agrees with H^ it is enough to recall that iV^ is defined as iVk restricted 
to its natural domain "K^ (]R^,P) fl Tir and to use the spectral calculus. These arguments 
justify the following: 



Proposition 3.11. The Harailtonian H , defined by (26), agrees on its domain of defi- 



nition with the Weyl quantization of the operator-valued symbol Hq defined by ([36j) . 
With a little abuse of notation, we refer to this result by writing H^ = Op^{Hq). 

3.6 Main result: effective dynamics for weak magnetic fields 

Let A^ and B^ be e-dependent (possibly unbounded) linear operators in "H. We write 
A,; = Bs + Oo{£°°) if: for any A^ G N there exist a positive constant Cm such that 

\\Ae-B,h^n)<CNe'' (40) 

for every e G [0,£:o)- Notice that, though the operators are unbounded, the difference is 
required to be a bounded operator. 

We refer to Appendices A and B of |Teu03| for the basic terminology concerning pseu- 
dodifferential operators, and in particular as for the notions of principal symbol, asymptotic 
expansion, resummation, Moyal product. 

Theorem 3.12. Let Assumptions (A^), (B), (C) and (D) be satisfied and let {-E'„(-)}nex 



(with \I\ = m) be an isolated family of energy bands for Hq satisfying condition (38). Then: 



1. Almost-invariant subspace: there exist an orthogonal projection He G ^{I-Lt), with 
IIe = O^i^iji) + (9o(^°°) (^'f^d the symbol Ti{k,vi) x Y1T=q^'^ '"'ji^^v) having principal part 
Tioik, rf) = Pxik — Airf)), so that 
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In particular for any A^ G N there exist a Cn such that 

||(l-n,)e-^^"n,||<C^e^|t| (41) 

for e sufficiently small, t G M. 



2. Effective dynamics: let "Href = L^{Mr*,cPk) ® "Hf, vTr as defined above (32) and 
Hi. = '^L'^[M-p*) ® ^r ^ ^i'Hrei)- Then there exist a unitary operator 

Ue '■ "Hr — >• "Href 

such that 

(i) Ue = Ope{u) + Oo{e°°), where the symbol u x YlT=o^^'^j ^'^■^ principal part Uq given 



by (32) with k replaced by K{k,ri); 

(ii) n, = Us n, u,-'; 

(iii) posing /C := Hi. "Href; one has 

u, n, H^ n, U-' = Hi^ + Oo(£~) G ^(/c) 

with if|g = Op^(/i) and /i a resummation of the formal symbol m jl tt jj iJo tl ttjjm"^ 
(^i/iits algorithmically computable at any finite order). Moreover, 

||(e-|^"_f/-ie-f^.^fff/^)n,|| <C^ 5^ (£+|t|). (42) 

Remark 3.13. The previous theorem and the following proof generalize straightforwardly 
to any dimension c? G N. We prefer to state it only in the case rf = 2 in view of the 
application to the QHE and of the comparison with the results in Section 3, the latter being 
valid only for d = 2. ♦(> 



Proof of Theorem 3.12 



Step 1. Almost-invariant subspace 

The proof of the existence of the super-adiabatic projection is very close to the proof of 
Proposition 1 of |PST03a| . so we only sketch the strategy and emphasize the main differences 
with respect to that proof. 

First of all, one constructs a formal symbol tt x Xl^o ^'^ ^j (^'^^ Moyal projection) such 
that: (i) TrjJTr x vr; (ii) vr"'' = vr; (iii) ifotl^ ^ ^tl-^o where x denotes the asymptotic 
equivalence of formal series. 

The symbol vr is constructed recursively at any order j G N starting from ttq and Hq. One 
firstly show the uniqueness of vr (see Lemma 2.3. in |PST03b] ). The uniqueness allows us 
to construct n locally, i.e. in a neighborhood of some point zq := (/cc^/o) G M^ x M^. From 
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the continuity of the map k \-^ Hper{k) and the condition ( p3T| it follows that there exists a 
neighborhood W^q of /cq such that for every k G W^o the set {£n{k)}nex can be enclosed by a 
positively-oriented circle E(fco) C C independent of k. Moreover it is possible to choose S(A;o) 
in such a way that: it is symmetric with respect to the real axis; dist{'E{kQ),a{H-per{k))) ^ 
^Cg for all k G Ukg', Radius (Zl(A;o)) ^ Cr is bounded by a constant Cr independent of ko; 
E(A;o - 7*) = ^(ko) for all 7* G T*. 



With the notation of Section 3.4 we have Hq = H^p o j^ with H^(k,T]) := Hper{k) + (f){r]). 
Let A{kQ,rio) := S(A;o) + 4'{Vo) denote the translation of the circle S(/co) by 0(?7o) ^ I^ and 
pose A := A o j^. From the smoothness of it follows that there exists a neighborhood 
Uzo C M^ X ]R2 of 2jj such that dist(A(zo), o-(-f^o(2;))) ^ |Cg for all z G W^^. Moreover A(zo) 
is symmetric with respect to the real axis, has radius bounded by Cr and is r*-periodic in 
the variable k = k — A{ri) (see Figure II]). 

We proceed by using the Riesz formula, namely by posing 

'7rj{z) := — (p d\ Rj{X,z) on U^o 

2^ Ja(zo) 

where Rj{X, •) denotes the j-th term in the Moyal resolvent R{X, ■) = YlT=o^''^ji^y ') (also 
known as the parametrix) , defined by the request that 



iHoi-)-\lv)^Ri\,-) = tnt, Ri\,-)UHoi-)-\tv) = tv on U, 



ZO- 



Each term Rj is computed by a recursive procedure starting from Ro{X, ■) := {Ho{-)—Xtv) ^, 
as illustrated in |GMS91] . Following |PST03al equations (30) and (31)] one obtains that 

Rj{\, z) = -Ro{\, z) Lj{X, z) (43) 

where Lj is the (j — l)-th order obstruction for Rq to be the Moyal resolvent, i.e. 

(/7o(-)-AMtl ($^£"i?,(A,-)) =l«, + £^L,(-) + 0(5^+^). (44) 

\n=0 / 

At the first order L\ = — | {Hq, -Ro}fc„; with {-, ■}k.ri the Poisson brackets. 
The technical (and crucial) part of the proof is to show that 

for all j G N, with v{k, 77) := (1 + \k\'^). By means of the recursive construction each Rj{X, ■) 
inherits the special r-equivariance from the principal symbol Ro{X, •) = ((i?^0oj^)(-)— Al©)"^. 
The special periodicity in n of the domain of integration A(-) which appears in the Riesz 
formula assures also the special r-equivariance of each ttj{-). 

Since ||(9>j)(2)||b ^ 27vCr sup;,gA(^o) \\9z{Rj){Kz)\\i, (b means either ^{Hi) or ^{Hf,V), 
a G N^ is a multiindex and d'^ := dk^dk^d'^fdm)i ^^ need only to prove that Rj{X,-) G 
S^.^ {^{'HfjV)) n Sl^.^ (ef^('Hf)) uniformly in A. This is the delicate point of the proof. 
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First of all, from the definition of A(2;o) it follows that 

||i?o(A,z)||,^(«,) = [dist(A,a(/7o(^)))]"' ^ Vc, 

uniformly in A. Let a G N'', with |cr| = 1. One observes that d'^Ro{X, z) = —Ro{X, z)N^{X, z) 
with N'^{\z) := d^Ho{z) Ro{X,z). From the relation 



z z z z z 

and an inductive argument, it follows the chain rule 

d:Ro = Ro Y.^p,...p,^,N!k..N^.^-^ 

where /^i, . . . , [3\a\ G N"^, |a| := ai + . . . + 04, ujp^...j3.. = ±1 is a suitable sign function and 
the sum runs over all the combinations of multiindices such that /3i + . . . + /3|o| = a with 
the convention N^ = 1. The chain rule implies that Rq G 5*^.^ (^(Hf)) provided that 

ll^ril.«(Wf) = \\9zHo RolsgCHf) < Ca uniformly in A. 

The latter condition is true since \\{d'^HQ){k,ri) i?o(A, A;,?7)||.5^(^j) ^ (g ° Jk)^^,!]), for a 
suitable g{k,ri), r*-periodic in k and bounded in 77; the latter claim can be checked as in 
Proposition 3.10[ 



Similarly to prove that Rq G S*^.^ {^{'Hf,V)) we need to show that 

\\Ro Ky^^nuv) = Witn, - A,) Ro Ky^^n,) < C^ v'i-) 

uniformly in A. Since A^" is bounded on "Hf it is sufficient to show that ||(l^j— Ae)i?o('^; -2)||^(-Hf) ^ 
C'^v'iz). Observe that ||(l«,-Ae)i?o(A, [K]-Y;v)\\.nn,) = \\{ln-^e)T{Y)-'Ro{X, M;r/)||,^(^,). 
The commutation relation 

-Aer(7*)-i = r(7*)-i {\^*\^ + i2^* ■ Ve - Ag) 

and the straightforward bound 

II (|7*P + ^2^* ■ We - A,) (1^, - Ae)-%j^n,) ^ 0(1 + I7I') < C'il + |«:(A;,r/)n 

imply 

11(1^,- Ae)i?o(A,^)||.^(«,) < C'^ v'izmtn,-Ae)RoiX,[K];r])yj^n,) 

with v' := V o j^. Finally observe that 

Uln, - Ae)Ro{X, [n];r])y(^n,) ^ fiH,v) ^ C" . (45) 

The first inequality above follows by an expansion on the Fourier basis, for fixed [k] and 
7]] the second follows from the fact that [k] takes values on a compact set and the explicit 



dependence on 77 is through the bounded function 0. The bound (45) implies that Rq G 
Sl.^ {^{ni,V)) n Sl.^ (=^CHf)) uniformly in A. 

To prove that Rj G S^.^ {^CHi}), we observe that for any a G N*^ one has 

d'^Rj{X,z) = Ro{X,z) M^.j{X,z) 
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where M^- is a linear combination of terms which are product oi N^ and d^Lj with |/3|, \6\ ^ 
|a|. Thus it is sufficient to prove that Lj G S^.^ (ef^(7/f)) for every j G N. The latter claim is 



proved by induction on j G N. Referring to (43), one has trivially that Li G S^.^ {^{'Hf)). 
Lj+i is a linear combination of products of N^ (with ^ |a| ^ j ' + 1) and M^- (with 
\P\ + i = j + 1 and ^ i ^ j). Then the induction hypothesis on Lj for alH = 1, . . . , j 
implies that L^+i is in 5*^.^ (=^('Hf)). 

Finally observing that ||(9"-Rj||.^(^j_x)) ^ ll^"jl|.«(Wf) ||-Ro||,:^(-Hf,i') ^ind using the fact that 
Ro G Si:.^ {^{Uu V)) it follows that R^ G Sl.^ (^(-Hf, V)) n S'^., (^{Ui)) uniformly in A, for 

all i en. 



As explained in Section 3.4, we can apply the result of Proposition B.4 in |Teu03) to spe- 
cial r-equivariant symbols obtaining Hq'^tt & S'^.^(^('Hf)). However the r-equivariance 
of ifo tl ^ and its derivatives implies that the norms are bounded in z, hence Hq'^'k G 
S\.^{li^{'Hi)) which implies by adjointness also Trjji/o ^ ^l.-ri'^i'^i))- ^Y construction 
[H^;0p^{7f)] = Op^{Ho^TT — TTJjifo) = C>o{e'^) where the remainder is bounded in the 
norm of ^{Tir)- 

The operator Op^{tt) is only approximately a projection, since Op^ln)^ = Op^lir^ir) = 
Opg(7r) + Oq{6°°). We obtain the super adiabatic projection He by using the trick in |NS04] . 
Indeed, one notices that, for e sufficiently small, the spectrum of Op^{tt) does not contain 
e.g. the points {1/2} and {3/2}. Thus, the formula 

n,= f / (Op,(7r)-.)-^ (46) 

yields an orthogonal projector such that Ilg = Op^{tt) + (9o(^°°)- 



Finally, equation (41) follows by observing that [H^; Ilg] = Oq^s'^) implies 

[e-i^";n,] = Oo(e°°|t|) 
as proved in Corollary 3.3 in |Teu03] . 

Step 2. Construction of the intertwining unitary 

The construction of the intertwining unitary follows as in the proof of Proposition 2 of 
|PST03a] . Firstly one constructs a formal symbol u x YLT=o^'''^3 such that: (i) v) '^u = 
u\\u'^ = Im; (ii) u\\n\\u'^ = n,. 

The existence of such a symbol follows from a recursive procedure starting from mq and 
using the expansion of vr x X^ylo ^''''^j obtained above. However, the symbol u which comes 
out of this procedure is not unique. 



Since Uq is right T-covariant (see the end of Section 3.3) in k, then one can prove by 



induction that the same is also true for all the symbols Uj and hence for the full symbol 
u. Finally, since uq G S^{^{'Hi)) one deduces by induction also Uj G S'^(=^('Hf)) for all 
j G N. The quantization of this symbol is an element of ^{l-Lr, ^ref) satisfying the following 
properties: 
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(i) Op,(w)0p,(w)t = l^^^, + Oo(£°°), 
(ii) Op,(w)tOp,(n) = l^^ + Oo(£°°), 

(iii) Op,(«)n,Op,Ht = n, + Oo(£°°). 

Nevertheless Oy>i,{u) can be modified by an OQ{e°°) term using the same technique of 
Lemma 3.3 (Step II) in |PST03b ] to obtain the true unitary Ug,. 

Step 3. Effective dynamics 

The last step of the proof is identical to the corresponding part (Proposition 3) of |PST03a] . 

3.7 Hofstadter-like Hamiltonians 

We now focus on the special case of a single isolated energy band £'*, i.e. m = 1, and we com- 
ment on the relation between the effective Hamiltonian, the celebrated Peierls substitution 
and Hofstadter-like Hamiltonians (see Section 1). 

In this special case, tto{k,) = | '?/'^, ( /t) )(■?/'* (k)! and Uo{k,) = \x){'ip*{'^)\ +'^0' where ip^:{k) is 
the eigenvector of Hperik) corresponding to the eigenvalue £^{k). Let h G S^{^{'Hi)) be a 
resummation of the formal symbol u'^'k'^Hq'^'k'^u^^. A straightforward computation yields 

ho = UQ'KQHQTiQul = \x){iljA |^*)(^*| -ffo |^*)(^*| |^*)(xl = -S*7r; 



* "r- 



Since tTj. is one-dimensional, /iq can be regarded as a scalar-valued symbol with explicit 
expression 

/io(A;, v) = E,{k, T]) = £,{k- A{r])) + <i){r]). 

By considering the quantization of the latter, the effective one-band Hamiltonian reads 

Op,(/io) = E,{k,ieVk) =£.{k- AiieVk)) + (t^iieVk)- (47) 

The latter formula corresponds to the momentum-space reformulation of the well-known 
Peierls substitution [Pei33l RM76] . 

To illustrate this point, we specialize to the case of a uniform external magnetic field and 



zero external electric field, setting = and Aq = in (36). The leading order contribution 
to the dynamics in the almost invariant subspace is therefore given by a bounded operator, 
acting on the reference Hilbert space L^(Mr*, rf^fc), defined as the quantization (in the sense 



of Section 3.5) of the function <£^* o j^ : (/c, rf) \-^ £^,{k — A{vi)), defined on T x 



d XX rod 



Loosely speaking, the above procedure corresponds to the following "substitution rule": 
one may think to quantize the smooth function £^, : T'' — >■ M by formally replacing the 
variables (/ci, /C2) with the operators (=^', =^') defined by 

Xi := k^ + ^(^,^)^, ^i := k^ - \{i,e) — , (48) 
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regarded as unbounded operators acting in L"^ {Mr* , d"^ It) . To make this procedure rigorous, 
one can expand £^* in its Fourier series, i.e. £^:{k) = ^nm&i^ri,m e*^'^(""+™^)''^' and define 
the Peierls ' quantization of S^, as the operator obtained by the same series expansion 
with the phases e'^'^(""+™^)'^ replaced by the unitary operators e*^'^(""+™''')''^' (the series is 
norm-convergent, in view of the regularity of £^*). This fixes uniquely the prescription for 
the quantization. 

To streamline the notation, one introduces new coordinates ^i := 27T{a-k) and ^2 := 27r(6- 
k) such that the function S'^^, ^^(1^1,^2) := <^*(^(0) becomes (27rZ) ^-periodic. The change 
of variables induces a unitary map from L^(Mr*^, (Pk) to L^(T^, (P^) which intertwines the 



operators (48) with the operators (recall e = "^'"/hb) 



so that 27r(a • JT') ^ Xi and 2Tx{h ■ J^f') ^ J^. 

Let F : T2 ^ C be sufficiently regular that its Fourier series F(^i, ^2) = E„,mez /n,»ne^("«i+'^«2) 
is uniformly-convergent. We define the Peierls quantization of F as 



I n,m 



e 



Let Uq = e*'^i and Vq = e*'^^ {Hofstadter unitaries), acting on L^(T^,d^^) as 

(lIo^)(ei,6) = e^^^^ Uu^2 - vr^") , (Vo^)(ei,6) = e^«^^ (^i + vr^,6 ) • (50) 



We regard (50) as the definition of the two unitaries, so there is no need to specify the 



domain of definition of the generators (49). Thus the Peierls quantization of the function F 



defines a bounded operator on L (T , d^) given, in terms of the Hofstadter unitaries, by 



+00 , 

Lq 



F{Uo,Vo)= Yl /n.^e^™"^"^^ ^o^^> (51) 



n,m=— 00 



where the fundamental commutation relation IXoVq = e ^ ^\'^bJ VqUq has been used. For- 



mula (51) defines a Hofstadter-like Hamiltonian with deformation parameter ^/hB- Indeed, 
the special case -^Hof = 1^0 + '^0^ + ^^^0 + ^0^ i^ (^P ^^ ^ unitary equivalence) the celebrated 
Hofstadter Hamiltonian |Hof 76| . 

Summarizing, we draw the following 
Conclusion 3.14. Under the assumption of Theorem ^.7, for every Vp G Lf^^{E?',d'^r), 



in the Hofstadter regime (Kb — ;■ ooj, the dynamics generated by the Hamiltonian Hbl ^ in 
the subspace related to a single isolated Bloch band, is approximated up to an error of order 
^/hs (and up to a unitary transform) by the dynamics generated on the reference Hilbert 
space L^(T^,(i^^) by a Hofstadter-like Hamiltonian, i.e. by a power series in the Hofstadter 



unitaries Uq and Vq defined by (50). 
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4 Space-adiabatic theory for the Harper regime 
4.1 Adiabatic parameter for strong magnetic fields 

We now consider the case of a strong external magnetic field. Since we are interested in the 
limit B — > +00 we set Aq = and $ = in the Hamiltonian dol). By exploiting the gauge 
freedom, we choose 

Vr- ■ Ar = 0, / Ar(r) £r = 0, (52) 

JMr 

this choice being always possible |Sob97| . Let us denote by Qr = {Qri,Qr2) the multipli- 
cation operators by ri and r2 and with P^ = {Pri,Pr2) = —ifiVr- Taking into account 
conditions (52) and Aq = 0, (p = 0, the Hamiltonian ^ is rewritten as 



H 



BL 



2m 



P, 



qB 



r\ T" r) Wr-2 



''2 r) Vri 



+ Vr(Q.)+W(Q, 



(53) 



where 



Vr(g.) = Vr(g.) + 



2mc^ 



|Ar(Q.) 



W(Q„P,) = -^(Ar)i(g. 



mc 



P +^Q ■ 



- — (Ar)2(Q. 



mc 



p — 



f2 






Vr-i 



(54) 
(55) 



with (Ar)i and (Ar)2 the F-periodic components of the vector potential Ap. The first of 



(52) assures that W is a symmetric operator. 

It is useful to define two new pairs of canonical dimensionless operators: 

1 ^^ 



(fast) 



K^:-- 



Ko:-- 



1 

25 
1 

25 



n 



(slow) 



a* ■ Qr ~ tg- b- P, 
n 



Gi :— - b ■ Qr -^ Lq^ Cl ■ Pr 

2 n 

G2 ■= - a* ■ Qr + Ln"-Tr b ■ P, 



i 

52 



(56) 



h 



where 5 := ^fhs = a/*o/Wz$^ according to the notation introduced in Section KLJ Since 
6"^ oc i/s, the limit of strong magnetic field corresponds to 5 — )■ 0. We consider 6 as the 
adiabatic parameter in the Harper regime. A direct computation shows that 



[K,; K2] = tLg In 



phy ' 



[Gi; G2] = It, 6' In,,,, [Kj- G^] = 0, j, A; = 1, 2. (57) 



These new variables are important for three reasons: 

(a) they make evident a separation of scales between the slow degrees of freedom related 
to the the dynamics induced by periodic potential and the fast degrees of freedom 
related to the cyclotron motion induced by the external magnetic field. Indeed, for 
Vr = 0, the fast variables {Ki,K2) (the kinetic momenta) describe the kinetic energy 
of the cyclotron motion, while the slow variables {Gi,G2) correspond semiclassically 
to the center of the cyclotron orbit and are conserved quantities. 
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(b) The new variables are diniensionless. According to the notation used in Section 3.1 
let ii/^BL '■= V^oHbl be the diniensionless Bloch-Landau Hamiltonian with £q := ^^/mnY 
the natural unit of energy. 

(c) The use of the new variables simplifies the expression of the F-periodic functions 
appearing in Hbl- Indeed, a* ■ Q,. = G2 + 5 K2 and h* ■ Qr = Gi — 5 Ki, hence if /r is 
any F-periodic function one has 

HQr) = f{G2 + 6K2,G,-6 K{) (58) 

where / is the Z^-periodic function related to /p. 



In terms of the new variables (56), the Hamiltonian H-bi, reads 



i^BL = ^ H(i^i, K2) + V{G2 + 6 K2, G,-6Ki) + ^ W{Ki, d, K2, G2) (59) 

where 

EiK,, K2) := ^ [\a\^K2' + \h\^K,' - a ■ b {K,; K2}] (60) 

is a quadratic function of the operators Ki and K2 ({■;■} denotes the anticommutator), 
V is the Z^-periodic function related to the F-periodic function Y^oVr and W denotes the 
function ^/eqW with respect the new canonical pairs, namely 

W{Ki, Gi, K2, G2) = h {G2 + 6 K2,G,~6 Ki) K^ - /s (G2 + <5 7^2, Gi - 5 K^) K2 (61) 

where /i and /2 are the Z^-periodic diniensionless functions 

/i(a*-r,6*-r) :=27r— ^(a*-Ar)(r) and ^(a* ■ r, 6* ■ r) := 27r— -^(6* ■ Ar)(r). 



An easy computation shows that the first gauge condition of (52) is equivalent to 



TT— (xi, X2) + 7^— (a;i, X2) = 0. (62) 

OXi 0x2 

Obviously Pi^ is a symmetric operator, since W is symmetric. 

The problem has a natural time-scale which is fixed by the cyclotron frequency ojc = -^^■ 



With respect to the (fast) ultramicroscopic time-scale r := oOcS, equation (12) becomes 



.i^^ = /.erf., *==!-. (63) 

Thus the physically relevant Hamiltonian is 

Hi^ := 6^ Hbl = 2(7^1,7^2) + S W{K,,G,,K2,G2) + 6^V{G2 + 6 K2,Gi - 6 K,) . (64) 
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4.2 Separation of scales: the von Neumann unitary 



The coniniutation relations (57) show that (i^i, K2) and (Gi, G2) are two pairs of canonical 
conjugate operators. The Stone-von Neumann uniqueness Theorem (see |BR97| Corollary 
5.2.15) assures the existence of a unitary map W (called von Neumann unitary) 

W : 7/phy — >n^ ■='H^®'Hi = L2(M, dx^) ® L2(M, dxi) (65) 

such that 

WGiW-i = Qs = multiplication by Xs, WGsW-^ = P^ = -^^7;^ (66) 



d_ 

dx 

The explicit costruction of the von Neumann unitary W is described in Appendix 



WKiW-^ = Qf = multiphcation by Xf, WfCsW"^ = Pf = -iia t^- (67) 

OXi 



Let Xj := Gj + {-iy6 Kj with j = 1, 2. From (66) and (67) it follows that 



X[ := WXiW-i = Qs-S Qu X'^ := WXaW'^ = Ps + S Pf. (68) 

Since Xi and X2 commute, one can use the spectral calculus to define any measurable func- 
tion of Xi and X2. For any / e L°°(M^ d^x) one defines /(Xi, X2) := ^2 f{xu X2) ^eIV ^2^ 
where dWj^^^ is the projection- valued measure corresponding to Xj. In view of the unitarity 
of W, and observing that dE''"'-^ := WdWi^^^W^^ is the projection-valued measure of X'-, one 
obtains that 

W/(Xi,X2)W-i= / f{xuX2) dE'J'UE'J'^ = f{X[,X',). 

So the effect of the conjugation through W on a function / of the operators Xi and X2 
formally amounts to replace the operators Xj with X'- inside /. 

In view of the above remark, one can easily rewrite H^^ making explicit the role of the 
fast and slow variables, obtaining 

H"^ := WHi^W-' = In^ ® E{Qu Pf) + 6 W{Qi, Q,, Pf, P,) + 6^V{P, + 6 Pf, Q, - 6 Qf) 

(69) 



where, according to (61), 



W{Qi, Qs, Pf, Ps) = h (Ps + 6PuQs-6 Qi) Qi - /2 (Ps + 5 Pf, Qs - 5 Qi) Pf- (70) 



4.3 Relevant part of the spectrum: the Landau bands 

The existence of a separation between fast and slow degrees of freedom and the decompo- 
sition of the physical Hilbert space 'Hphy into the product space "Hw = 'Hs^'Hi are the first 
two ingredients to develop the SAPT. According to the general scheme, we "replace" the 
canonical operators corresponding to the slow degrees of freedom with classical variables 
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which will be re-quantized "a posteriori". Mathematically, we show that the Hamiltonian 
H^ acting in "H^ is the Weyl quantization of the operator- valued function (symbol) H^, 

Hs{p.,x,):=E{Qi,Pf) + 6 W{Qux,,Pf,Ps)+6' V {p, + 6 P^x, - 6 Qf) . (71) 

V V 

The quantization is defined (formally) by the rules 

Xs I ^ Op5(Xs) := Qs ® 1%, Ps I ^ Op5(ps) := P^ ® l^r 



For every [p^^x.^ G M , equation (71) defines an unbounded operator Hs{ps,Xs) which acts 
in the Hilbert space ?^f. To make the quantization procedure rigorous, as explained in 
Appendix A of [ PST03b] . we need to consider Hs as function from M^ into some Banach 
space which is also a domain of self-adjointness for Hs{ps,Xs). We take care of this details 
in the Section l434l 

To complete the list of ingredients needed for the SAPT, we need to analize the spectrum 



of the principal part of the symbol (71) as {ps,Xs) varies in M^. The principal part of the 



symbol, denoted by Ho{ps,Xs), is given by (71) when 6 = 0, so it reads: 



HoiPs, Xs) := S(gf, Pf) = ^ [\a\'P,' + \b\'Q,' -a-b {Qf, P,}] . (72) 

Since the principal symbol is constant on the phase space, i.e. Hq{ps,Xs) = S for all 
{Ps, Xs) G M^, we are reduced to compute the spectrum of S. As well-known (see Remark 



4.1 below), the spectrum of S is pure point with cr(S) = {A„ := (n + 1/2) : n G N}. We 
refer to the eigenvalue A„ as the n-th Landau level. 

The spectrum of the symbol Hq consists of a collection of constant functions a^i : M^ — t- M, 
n &N, an{ps,Xs) = A„, which we call Landau bands. The band 0"„ is separated by the rest 
of the spectrum by a constant gap. In the gap condition (analogous to ([3l])) one can choose 
Cg = 1. Therefore, each finite family of contiguous Landau bands defines a relevant part of 
the spectrun appropriate to develop the SAPT. 

Remark 4.1 (The domain of self-adjointness). We describe explicitly the domain of self- 
adjointness of Hq{ps,Xs). Mimicking the standard theory of Landau levels, one introduces 
operators 

i i i 

a ■.= ^— [{ai + ia2)Pi - (61 + ib2)Qi] = ^[z-aPi-z-b Qf] (73) 

a^= ^^ [(ai - ia2)Pi - (61 - ib2)Qi] = ^ [za Pi - z, Q,] , (74) 

where Za := \{ai — ^2) and Zh '.= \{bi — 262 )• It is easy to check that 

aa"^ = S(gf, Pf) + i, ^1^,, a'^a = E{Q,, P,) - i, ^ln„ [a; a^] = i. In,- (75) 

Without loss of generality, we suppose that Lg = 1. Let tpo be the ground state defined 
by atpo = 0. A simple computation shows that ^o{x{) = C'e"'-^"*"-'^'' , where C > is a 
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Figure 2: Structure of the spectrum of Hq. The picture shows a "relevant part of the spectrum" consisting 
of two Landau bands of constant energy A* and A*+i. 



normalization constant, and a G M, /3 > are related to the geometry of the lattice F by 
a := '^■y2|a|2 and /3 := ^r/2|a|2. Since ipo is a fast decreasing smooth function, the vectors 
ipn '■= {n\)~2[a'y'"ipo, with n = 0,1, . . ., are well defined. From the algebraic relations (75) 
it follows straightforwardly that: (i) aipn = \/n'ilJn-i', (ii) the family of vectors {^„}„gN is 
an orthonormal basis for Hi called the generalized Hermite basis; (iii) Sipn = A^^n; (iv) the 
spectrum of S is pure point with o"(H) = {A„ : n G N}. 

Let C G Tif he the set of the finite linear combinations of the elements of the basis 
{?/'„}neN- The unbounded operators a, a^ and H are well defined on C and on this domain a^ 
acts as the adjoint of a and S is symmetric. Both a and a^ are closable and we will denote 
their closure by the same symbols. The operator S is essentially selfadjoint on the domain 
C (the deficiency indices are both zero) and so its domain of selfadjointness J-" := 2^(2) is 
the closure of C with respect to the graph norm HV'lll := ||V'|||^f + ||S'?/'||^j- The graph norm 
is equivalent to the more simple regularized norm 



T 



\-W\\Hi 



The domain J-" has the structure of an Hilbert space with hermitian structure provided by 



the regularized scalar product {ip; ip)jr := {Sip; 'Bip)'^^ 



♦0 
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4.4 Symbol class and asymptotic expansion 

In this section, we firstly identify tlie Banacli space in wfiicfi tfie symbol Hs defined by 



( 71 ) takes values and, secondarily, we explain in which sense H^ is a "seniiclassical symbol" 



in a suitable Hormander symbol class. The main results are contained in Proposition |4.2[ 
Readers who are not interested in technical details can jump directly to the next section. For 
the definitions of the Hormander classes S^{^{'Hi)) and S^{^{J^ ■,?{{)) we refer to Section 

Proposition 4.2. Assume that Assumption (Ag) holds true. Then for all (ps, Xg) G M^ the 
operator Hs{ps, Xs) is essentially self-adjoint on the dense domain C cHf consisting of finite 
linear combinations of generalized Hermite functions, and its domain of self-adjointness is 
the domain T on which the operator Hq = S is self-adjoint. Finally, Hs is in the Hormander 
class S\^{T; Hi)). 

In particular, Hs{ps,Xs) is a bounded operator from the Hilbert space J-" to the Hilbert 
space Hf for all {ps, Xg) G M^. The proof of the Proposition 4.2 follows from the Kato-Rellich 



Theorem showing that for any (pg, Xs) G M the operator Hs{ps, Xs) differs from Hq by a rel- 



atively bounded perturbation. The latter claim will be proved in Lemmas 4.3 and|4.4|below 



In view of Assumption (As), Vp G C{^(]R^,]R) so its Fourier series 

Vr(r) = ^ Wn,m e''™ "•-e^^.'" '>*■' 

converges uniformly and moreover X]nm=-oo |'^|"^|^|°'^ |'U^n,m| ^ C^ for all a G N^. 



Let V be the Z^-periodic function related to y^oVr, as in Section 4.1, In view of (58) 
one has 



+ 00 



W-^(Xi, X2)W-^ = V{Ps + 6Pf, Qs - SP{) = Y, ^"."^ e^2-(«^=+"^Q=) e*2-'5(«^f-™Qf) (75) 



n,m=— 00 



£0 

with Vn,m '■= V^oWn.m and where we used the fact that fast and slow variables commute 



and [Qs! -Ps] = S^[Q{',Pf]- The operator (76) can be seen as the Weyl quantization of the 
operator-valued symbol 

+00 
Vs{Ps,Xs):= Yl w„,^ e^2"('^P=+'"^=) e^2"''("-^f-™^f) (77) 

n,m=— 00 

with quantization rule 

Lemma 4.3. Let Assumption (As) hold true. Then V5 G S^{^{Hf)) n S^{BS{J^;'Hi)). In 
particular Vs{psiXs) is a hounded self-adjoint operator on Hf for all {ps,Xs) G M^. 
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Proof. It is sufficient to show that V5 E S^{^{l-Li)) since 



in view of IIH ^ 



2. Let a := (ai,a2) ^ N^, then 

+00 



|5p?C^^fe,Xs)|[^(^^,)^(27r)H Yl 



n,m-\ — 00 



. (2vr)H 



for all {ps, Xs) G M^, as a consequence of the unitarity of e*^'^''("^f mQf)_ r^^^^ self-adjointness 
follows by observing that {Vn,m} are the Fourier coefficients of a real function. ■ 



Assumption (As) implies that the F-periodic functions a* ■ Ar and b* ■ Ap are elements of 
^(M^, M). By the same arguments above, one proves that the operators fj{Ps + SP{, Qs — 



6Pf), j = 1,2, appearing in (70), are the Weyl quantization of the operator- valued functions 

+00 

f '■•'■'(» X ) •= '^ f^^'^ pi27r(nps+mxs) ^i2TrS{nPi-mQi) 



E /.'. 



n.m=— 00 



J = 1,2 



(79) 



according to (78). The coefficients ^^Q^fn,rn are the Fourier coefficients of a* ■ Ar if j = 1 



and of h* ■ Ap if j = 2. Thus, equation (70) shows that the operator W{Qi,Qs,Pi,Ps) 

(80) 



coincides with the Weyl quantization of the operator-valued symbol 

Ws{ps,Xs) := ff\ps,x^) Qi + ff\ps,Xs) Pi, 
defined, initially, on the dense domain C. 

Lemma 4.4. Let Assumption (A^) hold true. Then f'^p e S\^{nf)) n S^^iJ'i'Hf)), 
for j = 1,2. For all {ps,Xs) G M?, the bounded operators fs iPs^^s) 0,1^^ self-adjoint while 
Ws{ps,Xs) is symmetric on the dense domain C and infinitesimally bounded with respect to 
E. Finally Ws E S\^{J^;nf)). 



Proof. As in the first part of the proof of Lemma 4.3, one proves that fg E S^{SS{l-Lf)) fl 
S^{^{3^; 'Hi)) and its self-adjointness. The operator Ws{ps, Xs) is a linear combination of Qf 



and Pf, which are densely defined on £, multiplied by bounded operators. Using (62) one 
checks by a direct computation that Ws{ps, Xg) acts as a symmetric operator on £. Since Qf 
and Pf are infinitesimally bounded with respect to S, then the same holds true for Ws{ps, Xs), 
{ps,Xs) E M^. The last claim follows by observing that 



p{j) 



d-ifr>Xi)y^^.,n,) 



p{j) 



fn XiE 



■^-iii ^ II V ■^ 



■^^1 II II Ac« f 



5 Wmnc), 



with J = 1,2 and Xf = Qi or Pf. Since ||Xf H-i||.^(^j) ^ C and f^/^ E S^^CHi)), the claim 
is proved. ■ 



Lemmas 4.3 and 4.4, together with the fact that Hq = S is clearly in S^{^{J^;'Hf)) 



imply the last part of Proposition 4.2 
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4.5 Semiclassics: the C(^^)-approximated symbol 

In this section we consider the asymptotic expansion for the symbol Hs in the parameter 5. 



The Fourier expansion (77) for Vs and the similar expression for Ws-, namely 

n,m=—oo 

suggest a way to expand the symbol H^ in powers of 5. By inserting the expansion e 



i2-KSIn^rn 



^•^ "^! In,rn -I with /^^ := nPf — mQf, in (77) and (81) and by exchanging the order of 
the series one obtains the formal expansions 

+00 +00 

Vsips, Xs) ~ Yl ^^ ^i+2(Ps, Xs) Wsips, Xs) ~ Y ^^ ^i+i(Ps' ^s) (82) 

i=o j=o 



where 



?1,7?1= — OO 



n,m=— oo 



In view of (62), one easily shows that the operators Wj are (formally) symmetric. 

The justification of the formal expansions above requires some cautions: (i) we need to 
specify the domains of definitions of the unbounded operators In.m'' and consequently the 
domains of definitions of Vj and Wj] (ii) we need to justify the exchange of the order of the 



series in the equations (77) and (81). 



As for (i), one notices that 

X nzh — mZa 

For all (n, m) G Z^ the operators /„_„ are essentially self-adjoint on the invariant dense 
domain £ (their deficiency indices are both zero). The powers In,m'' are also well defined 
and essentially self-adjoint on C, as consequence of the Nelson Theorem (Theorem X.39 in 
|RS75| ) since the set {ipn}neN of the generalized Hermite functions is a total set of analytic 
vectors for every In,m (see Example 2 Section X.6 of |RS75j ). The domain of self-adjointness 
for In.m.'' is the closure of C with respect the corresponding graph norm. 



The operator Vj{ps, Xs) defined by equation (83) is an homogeneous polynomial of degree 
j — 2 in a and a^. It is symmetric (hence closable) and essentially self-adjoint on the invariant 
dense domain C. Analogously, the operators 

^(1) ^(2) 

M^ =1 ^' rf(^) Of+ f(2) Pfl =/ ^ \o a + Q — a^] a .^ ^^hfj^±^bjn^ 

■^"n,m ■ -^n,m [J n,m ^i ^^ J n,m ^ ij ^n,m \))n,m, " T^ i/n,m " J ; i/n,m, • i— 

(86) 
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which appear in the right-hand side of equation (84), are essentially self-adjoint on C since 



the set of the generalized Hermite functions provides a total set of analytic vectors. Thus 
we answered to point (i). 

Since the generalized Hermite functions are a total set of analytic vectors for all /„ ^ 
then the series ^ 



+00 (i27r(5)J 
3=0 



r- 



In,m'''^ converges in norm for every ip G C. From this observation, 



c(i) 



(2) 



the fact that the series of coefficients Vn^m, fn,m and fn,m are absolutely convergent and that 
Qf end P{ leave invariant the domain C, one argues that for all ■?/' G £ the double series 
which defines VslpsjXs)^^ and Ws{ps,Xs)ip are absolutely convergent, hence the order of the 
sums can be exchanged. Thus the series appearing on the r.h.s. of (82) agrees with Vs 



(respectively, with Ws) on the dense domain C By a density argument, the equality in (82) 
holds true on the full domain of definition of Vg (which is T-C^) and Ws respectively. 

In view of the above, we write the "semiclassical expansion" of the symbol Hs as: 



+00 



HsiPs, Xs) = H + ^ (5^ Hj{ps, Xs), Hj{ps, Xs) := Wj{ps, x^) + Vj{ps, 



Xo 



^7) 



with Vt = 0. 



Proposition 4.2 shows that the natural domain for the full symbol Hs{ps^ Xs) is the domain 



J-" of self-adjointness of S. However, if we want to truncate the series (87) at the j-th order 



we must be careful in the determination of the domain of definition of the single terms and 



to control the remainder. Every term in the expansion (87) is essentially self-adjoint on C 
However, the j-th order term Hj is the sum of two homogeneous polynomials in Qf and Pf 
(or equivalently in a and a^), Wj of degree j and Vj of degree j — 2. Since Wj = if Ar = 0, 
one obtains 

' J, if Ar 7^ 

J -2, ifAr = 



deg Hj 



where deg Hj means the degree of Hj as a polynomial in Q{ and Pf. If deg Hj > 2 then 
the operator Hj is not bounded by the principal symbol S, and in this sense it cannot be 
considered as a "small perturbation" in the sense of Kato. Moreover, some other problems 



appear (see Remark 4.8). In order to avoid these problems, we truncate the expansion (87) 
up to the polynomial term of degree 2, i.e. up to order 6"^ if Ar = and up to order 6'^ if 
Ar^O. 

Hereafter let \] be the indicator function of the periodic vector potential, defined as 

0, if Ar ^ 

1, if Ar = 0. 



2(1+W 



Let -ff^(ps, Xs) := S -F Y.j=i ^^Hj{p^, Xs), namely 



+00 



if](Ps,Xs)=S + 52 ^ 



_,j27r(nps+"iXs) 



n,m,=— 00 
+00 



tu, 



i2-n5In,r. 



1 



t2'KY8'ln. 



H^ip,, x,) = E + 6 J2 e^'^("^=+'"^=) [M° ^ + 6{z27rMl^ + t;„,„l^, 



(88) 
(89) 
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We call H\ the approximated symbol up to order (5^*-^^''-'. As a consequence of the Kato- 
Rellich theorem we have the following result: 

Proposition 4.5. Under Assumption (As) there exists a constant 6q such that for every 
5 < 5q and for every {ps, Xg) G M^ the operator H^iPs, Xs) (both for [\ = or 1) is self-adjoint 
on the domain T and bounded from below. Moreover H^ G S^{^{J-' ;!-[{)). 



Proof. As proved in Lemma A. 2 a and a^ are infinitesimally bounded with respect to H. 
This fact and Assumption (Ag), which assures the fast decay of the coefficients Vn,m and 
gn,m. (see ([86|), imply that the operators 

+ 00 +00 



Yl ^"'^ 



^i2w(nps+mxs) 



{In, + t27i5I„ 



Y^ ^^2.inp^+m.^) ^^0 ^ ^ 6Vn,^tn,] 



n,m=—oo 



n,m=—oo 



are infinitesimally bounded with respect to S and are elements of 5'^(e^(J-'; Hf)). 

The operators In,rn' and M^ ^ are only bounded (and not infinitesimally bounded) with 

respect to S. First of all it is easy to check that Ha" i^Wut ^ 3||H-?/'||^j for every ■?/' G J-", 
where a* means a or a^. Then, for every t/^ G J-", 

(-2 



\j-n,m yWh, ^ jCtn.ml 



la^^ll^, + ||{a; a^Yi^U + l|a^'^ll«,)' ^ 27^(^2 + m') 



''' ''m% 



(90) 

where we used the inequality (a + /3 + 7)^ ^ 3(a;^ + /3^ + 7^), the identity {a; a'''} = 2S and 
the bound |a„,mP ^ '^^/n{n'^ — m?) with d^ := max{|a|, |6|}. Assumption (Ag) assures that 



'Jn,m e*^'^'^"^="'''"^=)/„^m^5 which appears in (88), is bounded by 



the operator 5 2t: X]nm=- 

by a constant 5^C , with C oc Xln m=-oo ^n,m(^^ + ra^), and is in S^{^{J-'; T-Li)). The claim 

for Hg follows from the Kato-Rellich theorem fixing 6q := C~^. 

The claim for H^ follows in the same way proving an inequality of the type (90) for 



M, 



1 

n,m 



Cn,ma^ + Cn,ma^ + 23f?((i„,„)H + Lq^{dn,m) whcrC C„,, 



C^n,myn,m 



and d„ 



'^n,mgn,m- Obscrvc that the series of coefficients gn^m decays rapidly, then also the serie Cn^m 
and dn,m have a fast decay and in particular are bounded. This implies that in the inequality 



of type (90) we can find a global constant which does not depend on n and m. 



It is useful to have explicit expressions of the first terms Hj, in terms of a and a^ . 



From equations (87), (88) and (89), using the Fourier expansion of the derivatives of V and 



g := yV2 {zafi + 2b/2)) it is easy to check the following: 



- Case 1: Ar = - In this situation 

H} = E + 6^H2 + 6^H3 + 6^H^ 
with 



H2{ps,Xs) = V{ps,Xs) t-Hi 



Ha{Ps, Xs 



V2 
1 - 

4 . 



D,\\V) 2S + DliV) a^ + Dl{V) a^' 



(91) 
(92) 

(93) 
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where D^ is the differential operator defined by Dz 



y, d ^ _d_ 

'^ dxs dps 



and Dz is obtained by 



replacing Za and Zf, with Za and zi,. Since V is real, Dz{V) = Dz{V), which shows that if 3 is 
symmetric. The explicit expression of the second order differential operator \Dz\'^ := D^oDz 

- — (\a\^ — " ' "^^ .,.0 52 



\D., 



-2a -6 



dxsdps dpi 



(94) 



Case 2: Ar 7^ - In this situation 



H^ = E + 6H^ + 6^H2 



with 



Hi{ps,Xs) = g{ps,Xs) a + g{ps,Xs) a' 
H2iPs,x,) = Vtn,-V2Dzig)E 



1 



V2 



Dzig) a' + D^ig) a 



(95) 
(96) 



In the computation of (96) we used the first of the gauge conditions (|52|) which assures that 

1 



DM 



V2ni 



dxs \dxs dpsj 



\b\ 



2df2 

dps 



is a real function. From the definition of g, /i and /2 it follows that 

Zf 
g{a* ■r,b*-r)= vrv^— [(Ar)i - z(Ar)2] (r), 



(97) 



namely g is the dimensionless Z^-periodic function related to the F-periodic function (Ar)i — 
^(Ar)2j up to a multiplicative constant. 

4.6 Main result: effective dynamics for strong magnetic fields 

Preliminary estimates on the remainder 

The difference 5^^ := Hs — H^ is a self-adjoint element of S^{^{J-';T-L{)), which we call the 
remainder symbol. To develop the SAPT for the Harper regime we need to estimate the 
order of the remainder symbol. The next result shows essentially that 



HsiPs, X,) = Hjip,, X,) + 0{S'^^+'^), IT, HsiPs, Xs) TT, = TT, Hjip,, Xs) TTr + 0{S' 



2(^+l)+l^ 



where 



TTr := ^\'iPh){i^k 



(98) 
(99) 



j=i 



is the projection on the subspace spanned by the finite family of generalized Hermite func- 
tions {ipk.y^i- In other words, the error done by replacing the true symbol Hs with the 
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approximated symbol H\ (which has order 2(t] + 1) in 5) is of the same order of the ap- 
proximated symbol, so in this sense H\ is not a good approximation for H^. On the other 
side, what we need to develop the SAPT is to control the operator tTj. H^ tTi-, which is well 
approximated by TTr H^ VTr up to an error of order 2(t] + 1) + 1 in 5. 

Proposition 4.6. Let Assumption (As) hold true. Then 9i^ has order 0((5^(^+^)), i.e. 
there exist a constant C such that \\'^^{pstXs)\\.3S(v,Hi) ^ C^^*-^^^-* for all (ps,Xs) E M^. 
Moreover \\D\\ Tir\\m(Hi) = hr ^\\\ss{Ht) ^ C(52(^+i)+\ for all {ps,Xs) G M^ i.e. ^hJ tt^, tt^ 9^J 
and [Uv^; vr^] are SS{l-Li)-valued symbols of order 0{5'^'^^^^^^'^). 

Proof. (Case t] = 1) The explicit expression of the remainder symbol is 

+ 00 



^]{P,,X,)=5^ Y, ^n,me*'" 



(nps+rnxs) 



n,m=—oo 



A2TT5In 



- l^f + i2Ti5h 



1 



{i2txY5'L, 



and from (|100[) it follows that \\'^]{ps,Xs)\\o:g{v,Ht) ^ ^'^Yl°^=-oo \vn,m\J^n,m with 



(100) 



A, 



n,m ■= sup \{0} 



-^i2TTSI„^, 



-(ln, + i2TT6In,m + -{127:^6^1^^, 



E-^ip 



-Hi 



Hi 



(101) 



since HV'IIj^ •= ll'^'V'llwf &iid T = S^^l-Li. The operators In,m are essentially self-adjoint on C 
and we denote their closure with the same symbol. Since the operators In,m^ are positive, 
we can consider the resolvent operators Rn,m '■= {In,m^ + l^f)^^- Let suppose that 



Sn,m,\0 ) ■ 



-J,2-K&In,'. 



Im + t2TT6In,n. + -{i^TTyS'ln 



R 



n,m 



^ m, (102) 



^(Hf) 



for all n,m eTj, with sup5C((5) < +00. Then equation (101) would imply 



Noticing that /„_„^ = a„,„^ a^ + an,m^ a^ + 2|a„,mp S 



and observing that HS^-*^! 



k'H i^(^j) 



2^_ 



.^(-Ht) 



1 and lla"'' S ^ 



|.«(Wf) 



2a/2, one deduces from inequality (102) that 



+00 



\ml{p.,x,)y^T,,ni)^Ci[ Yl \vn,n.\{\n\ + \m\y]6'aS)^C26'aS) 



\n,m=— 00 



for suitable positive constants Ci and C2. 

It remains to prove the inequality ( 102[ ) and the estimate on C('^)- By spectral calculus 
one has that Cn,mi5) = supig^(^^_^) \Zsit)\ < sup^g^ \Zs{t)\ =: ((5) where 



J2TTSt 



2x2 



Z5(t) := An'S 



1 + t2n5t - -{2TiSt) 



{2n5ty + 477252 
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After some manipulations and the change of variable r := 2Ti5t one has that 



GAt) 



47rM4 



2715 



2 r^ + 4r2cos(r) -8rsin(r) -8cos(r) + 8 

^ ~A < ^3- 



Thus C(5)' = 47r4<54 sup^gK^^lr) ^ Ati^C^5\ hence ||9^KPs,Xs)||,^(i,,«,) ^ C5\ This con- 
eludes the first part of the proof. 

Since ||fHj7rr||^(«f) ^ YlT=i ll^llV'fc,)(^fcJIUcHf), then it is enough to show that for any 
Hermite vector ipk the inequality ||9iJ|?/'fc)(^fc|||,«(-Hf) ^ Ck5^ holds true. Observing that 
\\^]Wk){'4^k\\\.ss{Hi) = \\^]^k\\Hv one deduces 



lim(5 ||9^5|^fc)(^fc|||^(^j) = lim 



+00 



'^n.m.^ 



i27r{nps-\-mXs) 



'^ {i27iys^-^ 



V'A 



n,m=—oo 



0=3 



Wf 



+00 



^-vr^ 
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5^ |t;n,™| ||/n,JV^fc||w, < -n'C'WA^'Mn, = -irVik + ^V- tt^C =: Q 



n,m=— CO 



where C" := Z]n,m=-oo l"^",™! I«n,m|^ is finite in view of Assumption (As). 

This shows that for all 6 G [0, (5o) (for a suitable 6o > 0) the norm ||9^J|'?/'fc)('?/'fc|||.^{Wf) 
is bounded by Ck6^ and so it follows that ||9^j7rr||.^(^f) ^ mC6^ with C := maxi^...^m{C'fe,}- 



WUI 



m{Hi) 



m-rr 



<5"r||/^(Wf)- 



Finally \\n,Dll\\,aj(^ni) = IK^^ItTj. 

(Case t] = 0) The proof proceeds as in the previous case. Divide the remainder symbol 
in two terms 9^° = D^q + 9^5 where: 

+ 00 

n,m=— oo 

+ 00 



n,rn=—oo 



The control of 9^? is easy, indeed ||9^'|'||.<3?(x),-Hf) ^ 2||9^?||,<^(-^j) ^ 4C(52 where C := J2'^'^=_^ \vn,m\- 
Moreover (with the same technique used for the case \] = 1), one can check that for any 
Hermite vector ipk the function ti{6) := Jj- ||9^5(Ps, Xs)4'k\\'Ht is bounded by a constant Ck > 
in a suitable interval [0,(5o)- This assures that ||9^5 ^r||.«(Wf) is of order 0{6^). 

To control DIq we need to estimate S„^m := || (e*^"^-^"'™ — l>^j — i27i6In,m) M^,m'^~^\\.^{nt)- 
Let i?^ „ be the resolvent {In,m+i^nf)^^ ■ It is easy to check that ||(J„^m+^lwf)^nm'='~^ll.«('Hf) 
is bounded by a linear expression in \n\ and \m\. Indeed, as proved in Proposition 4.5, both 



M°„^ and M^ ,^ are bounded by E. By spectral calculus || (e*^'"^-^"'™ — 1^^ — i2Ti5In,m) R'n,m\\.gs('H{) 
is bounded by the maximum in r of the function F^ij-) := 47r^(5^ r -2rsmT-2cosT+2 ^ 

The last part follows observing that M^^ is a linear combinations of a and a^ and 
so they act splitting a Hermite vector ^/'fc as c^^ipk-i + dtm'^k+i where for a fixed A; the 
coefficients depend on fn,m- To conclude the proof it is sufficient to notice that to (5) : = 
_^|| ^gi27r<5/„,„ _ "i^^ _ i27i6Inm) "^fcHwf is bouuded by a constant C/j > in a suitable interval 
[0,<5o). ' ■ 
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Derivation of the adiabatically decoupled effective dynamics 



We recall that the Weyl quantization of the symbol Hs is the Hamiltonian (69), namely 
Op^(if5) = H^ . As for the approximated symbol H^^ we pose H^ := Op^(if^). Both if^ 
and H^ are bounded operators from L^(]R, dxs) ® J-' to Ti^ := L^(]R, dxs) ® "Hf. 

Theorem 4.7. Let Assumption (As) be satisfied. Let {an{-)}nex, with I = {n, . . . ,n + m — 
1}, be a family of Landau bands for S and let tt^. := Xlnex l'^"-)(^"l ^^ ^^^ spectral projector 
of Hq = S corresponding to the set {o'n{Ps,Xs)}nex- Then: 

1. Almost- invariant subspace: there exists an orthogonal projection U^ G ^(T-L^), with 
nj = Op5(7r) + Oo{5°°), 7r(ps,Xs) x Yl%o^^ T^jiPs,^^), and no{p^,x^) = tt, , such that 

[H\ Ii\] = Oo(5°°), [if^; nj] = Oo(52(^+^)+^). (103) 

2. Effective dynamics: letU, := tn,®T^T e ^CH^) andK, := Ran 11^ ~ L^{R,dXs)(S>C"^. 
Then there exists a unitary operator U^ G ^{T-i^) such that 

(i) Ul = Op^(m) + Oo{5°^), where the symbol u x Xl^o ^^''^3 ^^^ principal part uq = 1^^.; 

(ii) n^ = u's n! uf; 

(iii) Let h^ in S^{^{l-Li)) be a resummation of the formal symbol u^ti^H^ ^-k^u^ and define 
the effective Hamiltonian by H^^ := Opg{h^). Since [iJ^jj; Ilr] = 0, if^g is a bounded 
operator on K,. Then 

U\ n; H^ Ii\ Uf' = Hl^ + 00(52^^+')+') G e^(/C). (104) 

2'. Effective dynamics for a single Landau hand when Ap = 0; Consider a single 
Landau band a*(-) = A*, so that tCj. = \iIj^:){iIj^:\. Then, up to the order 6"^, one has that 

h!s = X*ln. + S' V(Ps, Qs) + <5' y F(Ps, Qs) + Oo {S') (105) 

where V{Ps,Qs) '■= Op^iV) is the Weyl quantization of the I?-periodic function V{j)s,Xs) 
related to the T-periodic potential \/ ^ , while y(Ps,Qs) := C)P5(|-Dzp(V^)) is the Weyl quanti- 



zation of the function \Dz\ {y){ps,Xs) defined through the differential operator (94) 



The derivation of the effective dynamics when Ar 7^ will be considered in Section 4.8 



Proof of Theorem 4.7 



Step 1. Almost-invariant subspace 



As explained in the first part of proof of the Theorem |3.12 one constructs a formal symbol 



TT (the Moyal projection) such that: (i) TrjJTr x tt; (ii) ti^ = tt; (iii) iL^jJTr x t^^H^. Such a 
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symbol vr x X]^o'^"'^j ^^ constructed recursively order by order starting from ttq 
H\ and it is unique (see Lemma 2.3. in |PST03b] ). The recursive relations are 



TTn •= 7r„ + 7r„ 



where the diagonal part is 7r„ := -n^Gni^r + (l^j — 7rr)G'n(l-Hf — vTr) with 



Cjn '■- 



'n-1 



'n-1 



Y.^'-AnY.^' 



TCi 



'n-1 

Y.^' 



vr,- 



0=0 / \j=0 / \i=0 

The off-diagonal part is defined by the implicit relation [Hq] tt^^] 



VTr and 



(106) 



(107) 



-F„ where 



n-1 \ /n-1 

j=0 / \j=0 



(108) 



The uniqueness allows us to construct vr only locally and this local construction is ex- 
plained in the second part of Lemma 2.3 in |PST03b] . In our case we can choose a {psiX^)- 
independent positively oriented complex circle A C C, symmetric with respect to the real 
axis, which encloses the family of (constant) spectral bands {cr„(-) = A„}„gx and such that 
dist(A, cr(ifo)) ^ \ (see Figurepl). For all A G A we construct recursively the Moyal resolvent 

(or parametrix) R\X; ■) := Yl'^o^''^ji^'^ ') '^^ ^s^ following the same technique explained 
during the proof of Theorem |3.12 The approximants of the symbol vr are related to the 

approximants of the Moyal resolvent by the usual Riesz formula i^j{z) := — §j^d\ RAX; z) 

where z := {PsiX^) E M^. Some care is required to show (iii) since, by construction, -f/^^tl^ 
takes values in SSifKi) while vr jj if^ takes values in ^(J-"). To solve this problem one can use 
the same argument proposed in Lemma 7 of |PST03a| . 

The technical and new part of the proof consist in showing that vr G S^l^lHf)) fl 
S^{^{HuF)). The Riesz formula implies \\{dyj){z)\\^ ^ 27r sup^^g^ ll^"^5(^; ^)llb for all 
a G N^ (b means either !^{1-Li) or ^{TifjJ^) and 9" := d'^^d^^) since A does not depend 
on z. Then we need only to show that i?J.(A; ■) G S^{^{Hi)) H S^{^{Hi,J^)). The choice 



of A assures \\Rq[X]z) 



.^{-Ht) 



XI 



-Hf) 



l.^CHt) 



^ 2. Moreover d'^Rl{X;z) = for all 



a 7^ and this implies that R^ G S^^^ifhLi)) uniformly in A. Since ||i?Q(A;z) 



|.«(-Ht,^) 



Al 



-Hf) 



\:nni) 



^ 00 one concludes that Wq E S^{ 



1{l-LuJ^)) uniformly in A. 



By means of equation (43), one has R- = —R^L- where L- is the j-t/i order obstruction for 
i?Q to be the Moyal resolvent. In view of this recursive relation, the proof of W- G S^{^{'Hf)) 
for all j G N is reduced to show that L) G 5'^(^('Hf)) for all j G N. 



The first order obstruction, computed by means of (44), is 



L\{X;z) = 6-'[{Hl{z) - Xln,)^Rl{X;z) - tn,]i = H,{z) rI{X;z) - -{H; i?i(A; z)}^,,,,. 

Since S and Rq do not depend on 2; G M^ it follows that L\ = HiRq. The operator Hi is 
linear in a and a^ (with all its derivative) if t] = or ifi = if t] = 1. In both cases Hi (with 
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its derivatives) is infinitesinially bounded with respect to S (see Lemma A. 2). This shows 
that L\ G S\^{Hi)) (but not in S\^{Hi,J^)) if \ = 0). 

We proceed by induction assuming that L- E S^{^{l-Li)) for all j ^ m G N. The (r7i + 1)- 



th order obstruction L^+i can be computed by means of equation (44) and the Moyal formula 
for the expansion of jj (see equation (A. 9) in |Teu03] ). After some manipulations, one gets 

^Ui(A; -) = 7^ E ^^^ {d:id;:H^R'^ (A; ^) (a-^^!) (A; -)• 

ai+02+r'+i=rn,+l 



Since HrRl G S^i^iUf)) uniformly in A (see Proposition |4^ then LJ^^^ G S^i^iUi)), and 
this concludes the inductive argument. 

Finally to prove i?5 G ^^(.^(T^f,;^)), observe that \\d'^R)\\,;g^nt,T) = \\^ R^ {d'^ L])y^nt) ^ 
Ca||Hi?J)||,^CHj) ^ +CX) for all j, a G N. 

Remark 4.8. It clearly emerges from the proof that the order ^^(h+i) jg ^j-^g ]-,gg^ g^p_ 
proximation which can be obtained with this technique. The obstruction is the condition 
HrRQ G S^(^{'H{)), which can be satisfied by the resolvent Rq := (S — C^f)^^ only for 
0^r^2(t] + l). 40 

Proposition A. 9 of |Teu03) assures that H^'\\n E S^{^{'Hf)) and, by adjointness, also 
vrJ^J G S\^{H{)). By construction [H^;Ops{'^)] = Op5([^J; 7r]j) = Oo{6°^) where 
[if^; 7r]jj := if ^ jj tt — tt jj if ^ = 0{6°°) denotes the Moyal commutator. Observing that 



[ii^; 7r]j = [W^ + 5^5;7r]j = [^s^ir]^ + C((5°°) and since Proposition |4.6| implies [9^^;7r]j = 
[9^5; T^j + 0(52(^+i)+i) = C(52(^+i)+i), it follows [ii^; 7r]j = 0(52(^+i)+i) which implies after 
the quantization [ii^; Op5(7r)] = Oo{6^'^^+^'>+^). 

The last step is to obtain the true projection 11^ (the super adiabatic projection) from 



estimates (103). 



Op^(7r) by means of the formula (46). Since 11^ — Opg(TT) = Oo{6°°), one recovers the 



Step 2. Construction of the intertwining unitary 

The construction of the intertwining unitary follows as in the proof of Theorem 3.1 of 
|PST03b] . Firstly one constructs a formal symbol u x Yl^o^'''^j such that: (i) u'^^u = 
u'^\u^ = 1-Hf ; (ii) m jj vr jj m^ = vTr. 

The existence of such a symbol follows from a recursive procedure starting from Uq 
(which can be fixed to be 1-^j in our specific case) and using the expansion of vr x Y1T=q ^'''^j 
obtained above. However, the symbol u which comes out of this procedure is not unique. 
The recursive relations are 

Un := a„ + hn with a^ := --^n, &n := K; Bn] (109) 
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where 



^n '■- 



'n-1 



'n-l 



5^5%, h 5^5^- 



Ua 



^-Hi 



0=0 



0=0 



and 



-Bn 



'n-l \ /n-1 

^ 5^Uj + 5"a„ U^ [Y^ ^'^J + ^" 

0=0 / \i=0 



vr^ 



;iio) 



:iir 



Since uq = 1^^ G S'^(=^('Hf)), then it follows by induction that Uj G S^i^CHi)) for all j G N. 

The quantization of u is an element of ^{l-i^) but it is not a true unitary. Nevertheless 
Op^(m) can be modified by an Oq{5^) term using the same technique of Lemma 3.3 (Step 
II) in |PST03b] to obtain the true unitary f/^. 



Step 3. Effective dynamics 

By construction [Hi^- H,] = Ov,{[h\ n,]^) = [U\ li] H\ hJ Ul^- H,] = since H, = U\ hJ U\ 



1-1 



Moreover equation (104) follows observing that U^ Yr^ H^ 11^ Ug 



H^Q coincides with the 



quantization of « jj tt jj 9^^ jj tt jj m"!" which is a symbol of order 0{5 



2(ll+l)+l^ 



Step 4. The case of a single Landau band when Ar = 

We need to expand the Moyal product h^^^ = -u jj vr jj H} tlvrtJM"l' = 7ri.tJ-utl H^ j| tt^ j] TTr + 0{5°°) 
up to the order 6^. To compute the various terms of the expansion /i^"^ x Yl'jLo ^^^j i^ i^ 
useful to define Xj •= ["^tl-^l tl'^^lj) so that hj = Ti^Xj'^r- Observing that 



'm—l 



m— 1 



n^H] 



Y.5^ xA ^u= iu^HlU^ -Y,^'xA U + 0{5n = 5'^Xm + 0{S 

0=0 / V i=0 



m+l^ 



one obtains the useful formula 



Xr. 



'm—l 



i=o 



'112) 



At the zeroth order one finds h^ = ttq Uq Hq Uq ttq = TTr S tTj. = A* n^ since Uq = l^j and 
ttq = TTi-. Its quantization is the operator Opg^ho) = A*l^^ acting on /C = L^(]R, (ixg). 

As for the first order (m = 1), Xi = ^i^o + UqHi - XoUi + K tl ^o]i - [Xo tl ^io]i = [^^i; -] 
since Xo = UqHqUq^ = H and Hi = 0. Then hi = 7ij.[ui; E]'Kj. = A*(7ri.ni7rr — tTj-UiTTj.) = 0, 
hence Opg{hi) = 0. 

At the second order (m = 2), one obtains after some manipulations X2 = H2 + U2'^ — 
'E.U2 — X1 "^1 which implies /12 = 'i^TH2'^T — '^TXi'^i'^r- We need to compute Ui. Using equations 
(109), (110) and (111) one obtains that — 2ai := \uq jJMq^ — l-Hfji = and hi := [tTi-; Bi\ with 
-Bi = [mo tl '^ tl "^0* — TTrji = TTi siucc tti = 0. To computc TTi wc use equations ( 106 ), ( 107) and 
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(108). Since Gi = [vTr U tTi- — 7rr]i = it follows that Trf = 0. In the case of a single energy 
band in the relevant part of the spectrum, the implicit relation which defines tt^^ can be 
solved, obtaining the useful equation 



vr" 



o° -= 7r,F„(S - Ktn,r\tn, - n,) - (In, - n,){E - Ktn^y'F^^r- (113) 



-Ki 


n,)-\lH, - n,) - 


H}h 


= Hi-K, 


- TT.iJi = 


r Ui = 


hi = 0. 


Then /la 



Since Fi = [Hi jJtTi. — 7ij.^Hg]i = HiTTj. — tTt-Hi = 0, being Hi = 0, it follows Bi = 7ii = 
TT^^ = and consecutively Ui = bi = 0. Then /i2 = TTj.H2TTr = VtTj., according to (91), and 
its quantization defines on /C the operator Opg{h2) = V{Ps,Xs). 



Considering (112) at the third order (m=3) and using Ui = 0, one obtains after some 
computations xs = H3 + u^E — Eu3 — X1U2 which implies h^ = tIj-HsTTj- — nT-XiU27ir. 
Thus we need to compute U2- Since ui = 0, it follows —2a2 = [uq jJuq^ — tnth = 0, B2 = 

[uq jj TT tj Uo^ — 7rr]2 = 7T2 and 62 = [TTr; TT2]- SiuCC TTi = 0, OUC haS that G2 = [vTr jj TTr — 7rr]2 = 

which implies tt^ = 0. To compute vr^^ we need F2 = [Hg^iij. — 7T^^HI]2 = [H2;tTj.] = 
[1-Hv ^r] = 0, where H2 = V l^j has been used. Then B2 = 112 = ^r^^ = and consequently 



U2 = 1)2 = 0. Therefore h^, = tTj-H^tt^, and equation (92) implies that tTj-H^tTj. = in view of 



TTj-aiTT- = {i/j^:\a\ip*)7Tj. = and similarly for a^ . Then Op^^h^) = 0. 

To compute the fourth order, we do not need to compute U3 and tt^. Indeed, by computing 



( 112 ) at the fourth order (m=4) one finds Xa = H^ + u^E — Eu^ + u^ Hi ~ xs^i = H4 + 
M4S — S'U4 since Hi = Ui = U2 = 0. Then h^ = ti^H^t^t- = ^\Dz(y)\'^ TTr, according to 



equation (93), and its quantization yields Opg(h4) = -fY{Ps,Q: 



4.7 Harper- like Hamiltonians 



The first term in (105) is a multiple of the identity, and therefore does not contribute to 
the dynamics as far as the expectation values of the observables are concerned. The leading 
term, providing a non-trivial contribution to the dynamics at the original microscopic time 
scale s oc 6'^t, is the bounded operator V{Ps, Qs) acting on the reference Hilbert space 
L^(R, (ixs). This operator is the Weyl quantization of the Z^-periodic smooth function V 
defined on the classical phase space M^. Hereafter we write Xg = x to simplify the notation. 

The quantization procedure can be reformulated by introducing the unitary operators 
Uoo := e"*^'^'^" and V^o '■= e~*^'^^= {Harper unitaries), acting on L^(R, (ix) as (recall 6"^ = 

'^b/2tt) 

(UooV^)(x) = e-*2^Xx), (V„oV^)(x) = ^{X - LghB). (114) 

For any Z^-periodic function F{p,x) = ^„m=-oo /"."iG"*^'^^"^^'"^'' whose Fourier series 
is uniformly convergent, the /i^-Weyl quantization of F is given by the formula 

+00 

n,m=— 00 

where the fundamental commutation relation Uno Voo = q-^'^'^^'-i^b) ^ H has been used. 



Formula (115) defines a Harper-like Hamiltonian with deformation parameter hs- Indeed, 
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the special case -ffnar = Uoo +^00^ + Voo + Vq^^ is the celebrated Harper Hamiltonian, namely 
the operator acting on L^(]R, dx) as 



In analogy with Section 3, we summarize the discussion in the following conclusion 



Conclusion 4.9. Under the assumptions of Theorem J^.l, for every Vr G C{f (M , M), in 
the Earner regime (Kb — ?• Oj, the dynamics generated by the Hamiltonian Hbl ^ restricted 
to the spectral subspace corresponding to a Landau level A=k is approximated up to an error of 
order Kb (and up to a unitary transform and an energy rescaling) by the dynamics generated 
on the reference Hilbert space L^(R, dx) by a Harper-like Hamiltonian, i.e. by a power series 



in the Harper unitaries Uqo and V^, defined by (114). 



4.8 Coupling of Landau bands in a periodic magnetic potential 



According to Theorem 4.7, the first non-trivial term which describes the effective dynamics 
in the almost invariant subspace related to a single Landau band A* is of order 6'^ oc Hb- 
An important ingredient in the proof is that Ar = implies Hi = 0. Moreover, the second 



non-trivial correction appears at order 6 oc Hb although H3 7^ 0. Indeed, the correction 



at order 6 vanishes since i^s, defined by (92), is linear in a and a\ hence {i/j^lH^li/j^) = 0. 



This observation suggests that for a family of Landau bands which contains two contiguous 

3 

bands {A^,,A*+i} one has, in general, a second non-trivial correction of order S^ oc Hb^ 
for the effective dynamics. Indeed, in this case one has irj-H^TTj. 7^ since {ijj^,\H3\ijj^^i) is 
generally non zero. Nevertheless, also in this case, the first non-trivial correction is of order 

Is there any mechanism to produce a non-trivial correction in the effective dynamics with 
leading order 5 oc /i^^ ? An affirmative answer requires Hi 7^ 0, and the latter condition is 
satisfied if we include in the Hamiltonian Hbl the effect of a F-periodic vector potential Ar 
(i.e. t] = 0). Since in this situation Hi is linear in a and a^, to obtain a non-trivial effect we 
need to consider a spectral subspace which contains at least two contiguous Landau bands. 

Our goal is to derive the (non-trivial) leading order for the effective Hamiltonian in this 



framework. According to the notation of Theorem 4.7, we need to expand the Moyal product 

f^\]=o _ ,y jj ^ jj jjo [JTrjj^it = vTi-jj-utJ H^ tj m"!" tj TTj. + 0{6°°) up to the first order 6. The symbols 
TT = TTr + (9(5) and u = l^j -|- 0{6) are derived as in the general construction showed in the 
proof of Theorem I^T) Now /C := Ran H^ ~ L'^{R, dx^) ® <C? . 

Expanding at zero order one finds Hq = -KqUqHqUq^ txq = tTi- H VTr = VTr S = S VTr and its 
quantization is the operator on /C defined by 

op,(ft„) ^ ( ("• \i) »- („, ;,, ,^^ ) ^ in. + 1)1. + ( i;- _,«^^ ) (117) 



As for the next order, from equation (112) it follows Xi = Hi + UiHq — HqUi (we use 



Xo = Hq) which imphes hi = 'n-rXi'^r = t^t Hi vTr + tc^ [ui; S] n^. To conclude the computation 
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we need Ml and vti. Using the recursive formulas (106), (107), (108), (109), (110) and (111), 



one obtains — 2ai := [mqIIwo^ — l^fj^ = 0, &i := [tTi-; i?i] and Bi = [mq tl ^ tl ""o* ~ '^t]i = tti 
since ai = 0. Observing that Gi = [TTrtJTTi. — 7rr]i = 0, it follows that vrf = and so 
Ml = [7rr;7ri] = [tTi-; tt^^] which implies TTi-UiTTr = 0. Finally 7rr[ui;S]7rr = 7rr[ui; TTrHTTrjTTi. = 



and so hi = tTj-Hi-k^. According to (95) the quantization of hi is an operator which acts on 
/C as 

/ s(Ps,gs)\ 

Opsihi) = V^^7+l (118) 

\9{Ps,QsV / 

where the operator S{Ps,Qs) is defined on L'^{M.,dXs) as the Weyl quantization of the Z^- 



periodic function g defined by equation (97). Summarizing, we obtained the following result: 



Theorem 4.10 (Effective Hamiltonian with a periodic magnetic potential). Under the 
assumptions of Theorem \4.1\ in the case Ap 7^ the dynamics in the spectral suhspace related 

,(■) = A*+j I j = 0, 1} is approximated by 

L^(R,dx,) (g)C2 



to a family of two contiguous Landau bands {a 
the effective Hamiltonian H^g := Op^{h^^ 
which is given, up to errors of order 5"^ , by 

1 



*^^) on the reference space K, 

S9{Ps,Q 



h: 



eff 



[n 



* + l)l/c + Vn^ + l 




1 



2Vn.+l 



1 



+ Oo {6') 



(119) 



Hs 



according to the notation introduced in (117[) and (118). 



Equation (97) shows that g{ps,Xs) = gi{Ps,Xs) — ig2iPs,Xs) where the function gi and 



g2 are related to the component (Ar)i and (Ar)2 of the F-periodic vector potential by the 
relation gj{a* ■r,b* -r) = 7rA/2|^(Ar)j(r), j = 1, 2. Let Sj(-Ps, Qs) be the Weyl quantization 
of gj. By introducing the Pauli matrices 



CTl 



1 

1 



0-2 







cr± 



(120) 



one can rewrite the effective Hamiltonian (119) in the form 

2 



H. 



eff 



(n, + l)lc2 + ^ax) ®l«. + (5vVTI ^(T,®gj(Ps,gs) + Oo('5'). (121) 
^ / i=i 



Clearly, the operator 9j{Ps, Qs) are Harper-like Hamiltonians and can be represented as a 
power series of the Harper unitaries IXoo and V^o of type ( 115 ). In this case the coefficients in 



the expansion are (up to a multiplicative constant) the Fourier coefficients of the components 
(Ar)j of the F-periodic vector potential. 

The determination of the spectrum of H^g can be reduced to the (generally simpler) 
problem of the computation of the spectrum of SS^- 

Proposition 4.11. Let Ht^^ be the first order approximation of the effective Hamiltonian 



'eff 



( |119| ) (or ^2\)). Then 

a{Hlf) = {n, + l) + S+VJS^, 



^±:={±v/V4 + (52(n, + l) A : A G a(ggt)} 
where a(9S'^) = a(ggt) u {0} if {0} G a(gtg) \ (7(ggt) and a(ggt) = (T(ggt) otherwise. 
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Proof. We give only a sketch of the proof. The term (n^, + l)lyc shifts the spectrum of a 
constant value {n^ + 1), then we can consider only the spectrum of 23 := H^^^ — {n^ + l)!^:- 
A simple computation shows that 

^-[6V^i;+i9' -11^, ) -ii^ + n^* + i)^ gtg 

which implies that cr(S2) = {1/4 + S^{n^ + 1)A : A G (t{S5^) U cx{S^Q)}. The operators 
SS^ S^S and S are bounded and self-adjoint. To show that (t(S9^) \ {0} = cr(S"''S) \ 
{0}, let A G cr(SS"'') with A 7^ and {ipn}neN C "Hg \ Ker(S^) be a sequence of non zero 
vectors such that ||(SS^ — ^)4'n\\Hs ~^ {Weyl's criterion), then ||(S^S — A)SVn||ws ^ 
||S^IU(w.)||(Sgt - \)i^n\\n. -^ 0. This implies that a(ggt) U a(gtg) = a(ggt). Now let 
e±{X) := ±a/V4 + 6'^{n^, + 1) A with A G a(gg"'') and {"ipnjnm a sequence of generalized 
eigenvectors for gg'l' relative to A. Then '^n ■= ((V2 + £±)4'n, ^a/w* + IQ'^ipn) G Hg ® C^ is 
a sequence of generalized eigenvectors for S relative to e-i-. ■ 



A Some technical results 

Since we include in our analysis a periodic vector potential Ap (as a new ingredient with 
respect to the standard literature) we include a short discussion of the self-adjointness and 
the spectral properties of the operators Hbl and if per- 



A.l Self-adjointness of Hbl and H, 



per 



The second Sobolev space "K (M. ) is defined to be the set of all ip G L (M ) such that 
d'^^d^^-ip G L^(]R^) in the sense of distributions for all n := {ni,n2) G N^ with 



\n\ 



ni + n2 ^ 2. One can proves that J{^(]R^) is the closure of C^(]R^,C) with respect to 
the Sobolev norm \\ ■ \\r^2 := ||(1 — A^-) ■ \\l^ and has a Hilbert space structure. Similarly 
the second magnetic- Sobolev space CK^(]R^) is defined to be the set of all tp G L^(]R^) such 
that Dl'^D^'^ii G L^(M^) in the sense of distributions for all n G N^ with |a| ^ 2, where 
Di := {dx^ + 1x2) and D2 '■= {8x2 — \xi). One can prove that IK^(]R^) is the closure of 
C^(M^,C) with respect to the magnetic- Sobolev norm || ■ ||j£2 := ||(1 — Aa/) • ||i2, where 
A^f := Di^ -\-D2^ is the magnetic- Laplacian. Moreover, !K^(]R^) has a natural Hilbert space 
structure. For further details see Section IX. 6 and IX. 7 of |RS75] and Chapter 7 of jLLOl] . 



Proof of Proposition 2.1 



We prove the claim for the dimensionless operators, namely we fix all the physical constants 
equal to 1 in (|9| and (15). 



- Step 1. First of all we prove that if per is essentially self-adjoint on C^(]R , C) and self- 
adjoint on D^2(]R2). Notice that 

^per = \ [-i^x - Ar{x)f + Vv{x) = -^A, + Ti + ^T2 (122) 
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with Ti := iAy ■ Vx- and T2 := ^(Va: ■ Ay) + |^rP + 2Vr- The free Hamiltonian — y2A^ is 
a self-adjoint operator with domain 3-C^(M^), essentially self-adjoint on C^(]R^, C) and from 
Assumption (A^) it follows that T2 is an infinitesimally bounded with respect to — y2Ax 
(notice that T2 — 2Vy is bounded). The symmetric operator Ti is unbounded with domain 
V{Ti) D J{2(M2). Let i) G 'K'^{M?), then 

with ip{X) the Fourier transform o{ip{x). For every a > 0, if 6 = ^ then A^ ^ (a|Ap + 6)^. It 
follows that ||(Ar)j9x-'?/'||j^2 ^ C'||(a|Ap + 6) 'i/'||^2 which implies that for all a' > (arbitrary 



small) there exists a 6' (depending on a') such that 

||(Ar)A.^lll. ^ « II 1^1' ^lli^ + ^' ll^lli^ = « IIA.^IIi^ + &' ll^lli- 

This inequality implies that Ti is infinitesimally bounded with respect to — y2Aa; and the 
thesis follows form the Kato-Rellich Theorem (see |RS75] Theorem X.12). 
- Step 2. The Bloch-Landau Hamiltonian is 

^BL := I HV, - Ar (x) - A {x)f + Vr (x) + <f){x). (123) 

Assumptions (A„) and (B) imply that {Ar + A)j e C^M^ M), j = 1, 2, and Vr + G Ll^CM."^) 
and this assures that Hbl is essentially self-adjoint on C^(M^, C) (see |RS75] Theorem X.34). 
Let A = Aq + Ab be the decomposition of the external vector potential with Aq smooth and 
bounded and Ab = ^{'-X2,Xi). By posing D := V^ — iAb, ^m '■= l-Dp and A^ := A^ + Aq, 
the Hamiltonian H^i^ reads 

Hbl = -^Am-A^-D + ^T. 

where T := i{'Vx-A\y) + \A\^\'^+2{Vr+4>). The operator T— 2Vr is bounded and the observation 
that y4b ■ -D is infinitesimally bounded with respect to — y2Aj\/ is an immediate consequence 
the assumption fj^^ \Vr{x)\'^ dPx < +cxd implies that Vr is uniformly locally 



A.2 



of Lemma 

L^ and hence, infinitesimally bounded with respect to — A^^ (see Theorem XIII. 96 in [ RS78| ). 
As proved in |AHS78j (Theorem 2.4) this is enough to claim that Vr is aslo infinitesimally 
bounded with respect to — A^/. Therefore, by the Kato-Rellich Theorem it follows that 
the domain of self-adjointness of i^BL coincides with the domain of self-adjointness of the 
magnetic-Laplacian, which is J{^(]R^). ■ 



A.2 Band spectrum of iifper 

We describe the spectral properties of the periodic Hamiltonian. The Bloch-Floquet trans- 
form maps unitarily H^er in H^^^ := J^^, Hper{k) <Ph- Then to have information about the 
spectrum of -ffper we need to study the spectra of the family of Hamiltonians 

Hver{k) = ^ [-tVe + k- Ar{e)f + Vr{e) = -^Ae + Ti{k) + ^T2{k) 
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where Ti{k) := i(Ar -\)-Ve and T2{k) := i(Ve ■ At) + \k\^ 
acting on the Hilbert space "Hf := L'^{Y, cPO) with V := M? /T 



+ l^rP + 2Vr are operators 
[Voronoi torus). 



Proof of Proposition 3.3 

- (i) The operator —^/2/\e on the Hilbert space Hf, is essentially selfadjoint on C°°(V), has 
domain of self-adjointness V := CK^(V) and its spectrum is pure point with {e*^''^*}^.gr* 
a complete orthogonal system of eigenvectors. If Assumption (A^) holds true then T2{k) 
infinitesimally bounded with respect to — y2A5), indeed T — 2Vr is bounded and Vr is in- 
finitesimally bounded (see |RS78] Theorem XIII. 97). With a Fourier estimate similar to 



those in the proof of Proposition 2.1 one can also show that Ti(k) is infinitesimally bounded 
with respect to — y2Ae, hence the Kato-Rellich Theorem implies that Hp^rik) is essentially 
self-adjoint on C°°(T^) and self-adjoint on the domain V. Moreover since — y2A6) is bounded 
below then also Hper{k) is bounded below. 

- (ii) For all ( in the resolvent set of —y2Ag the resolvent operator ro(C) := (""VsAg — (^l^J^^ 
is a compact operator. Since Ti{k) + ^T2{k) is a bounded perturbation of —'^/2Ag it follows 
that if per (/c) has compact resolvent (see |RS78j Theorem XIII. 68) and moreover it has a 
purely discrete spectrum with eigenvalues £n{k) — )■ +oo as n — )■ +oo (see |RS78| Theorem 
XIII.64). 

- (iii) The continuity of the function Sn{-) follows from the perturbation theory of dis- 
crete spectrum (see |RS78) Theorem XII. 13). Indeed, as discussed in Remark 3.2, ifper(-) 
is an analytic family (of type A) 



m 



the 



sense 



of Kato. Finally, since Hpei{k — 7* 



r(7*)ifper(fc)r(7*) ^, with r(7*) a unitary operator, then Sn{-) are r*-periodic. 



A. 3 The Landau Hamiltonian Hj^ 



The Landau Hamiltonian is the operator 



H, 



-^Am 



2 V 1 ^ 2; 2 



d 1 

oxi 2 



+ 



—i- 



d_ 
dx-. 



1 



(124) 



where Kj := —iDj, with j = 1,2, are the kinetic momenta. The Landau Hamiltonian i^L 
is essentially self-adjoint on on C~(]R2.c) c ^^(IR^) (gee |RS75j Theorem X.34) and its 
domain of self-adjointness is exactly the second magnetic-Sobolev space J{^(]R^) defined 



in Section A.l To describe the spectrum of H^ is helpful to introduce another pair of 



operators: Gi 



-id. 



XI 



^X2 and G2 '■= idx2 — hxi. The operators Ki,K2,Gi,G2 are all 



essentially self-adjoint on C^(]R^;C) (they have deficiency indices equal to zero) and on 
this domain the following commutation relations hold true 



[K,■,K2] = [G^■G2]=^1, [G,;Ki_ 



0. 



(125) 



The last of the (124) implies [Gj;Hi2 = hence the operators Gi and G2 are cause for 
the degeneration of the spectral eigenspaces of H^. It is a common lore to introduce the 
annihilation operator a := ^/V^{K2 — iKi) (its adjoint a^ is called creation operator) and the 
degeneration operator q := y'v^(G2 — iGi)- They fulfill the following (bosonic) commutation 
relation 



[a; a"^] = [0; 0"^ 



1, 



[0; Ht] = [0^; Ht] = 0, [a; H^] = a, [a^; H^ 



-at. 



(126) 
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The last two relations follow from the equality Hi^ = 00^ — 1/2! = 0^0+1/2!. Define the ground 
state ipQ G L'^iM?) as the normalized solution of QipQ = = aipo, i-G- V'ol^;) = Ce^jl^l . The 
generalized Hermite function of order (n,m) is defined to be 4'n,m '■= A — , (fl^) (a^) V'o- 
We will denote by J-" C L^(]R^) the set of the finite linear combinations of the vectors ipn^m 
and we will call it the Hermite domain. Clearly J-" C ^(IR^) (the Schwartz space). 

Lemma A.l. With the notation above: 

(i) the set {'ipn^m '■ n,m = 0,1,2, ...} is a complete orthonormal basis for L^(]R^) and so T 
is a dense domain; 

(ii) the spectrum of Hi is pure point and is given by {A^ := {n + 1/2) : ra = 0, 1, 2, . . .}, 
moreover Hi ipn,m = A„ ipn,m for every m = 0,1,2, .. . (degeneration index); 

(iii) Hl is essentially self-adjoint on J-" and the closure of T with respect to the magnetic- 
Sobolev norm coincides with the magnetic- Sobolev space !K^(]R^). 

Proof. - (i) Let W : L'^{M.^, d'^x) -^ L'^iR, du)'S)L'^(R, dv) be the unitary map such that the 
conjugate pairs {Ki, K2) and {Gi, G2) are transformed by W . . . W~^ into the canonical pairs 
(u, —idu) and {v, —idy). The existence of such a unitary W will be discussed in Appendix 
[BJ Obviously a ^-)• a = V\/2(m + du), i-^ = yV2{v + d^) and -^o := ^^ipo is the solution 
oi Qipo = aipQ = 0, namely ipo{u,v) = Hq^u) /io(^) where hoit) := 7r~4e~2* is the 0-th 
Hermite function. Then il)n,m{u,v) := (W-?/'„„^)(M,f) = hn{u) ® hm{v) which shows that 
the functions il'n,m define an orthonormal basis for L^(]R, du) ® -L^(]R, dv) since the Hermite 
functions /i„ are an orthonormal system for L^(]R). The claim follows since W is a unitary 
map. 

- (ii) Clearly Hi^tpo = {a^a + ^/2l)-ipo = ^/2ipQ and from relations ( |126 ) it follows that 

HLlpn,m = ;7= (fl"^)™^! (a"^)"^0 = V#n,m + ^7= (fl^)™ (o^Cl) (a'^)"^0 = A„ Vn,m- Then 

the generalized Hermite functions ^pn,m are a complete set of orthonormal eigenvectors for 
Hi^. This proves that the spectrum of Hi^ is pure point. 

- (iii) The operator Hi^ is essentially self-adjoint in J-" since the deficiency indices are both 
zero. This implies that last part of the claim. ■ 

Lemma A. 2. The operators Ki, K2, a and a^ are infinitesimally bounded with respect to 

Proof. Since Ki = 1/^2(0 + a') and K2 = ^/iV2{a — a') it is enough to prove the claim for 

a and a^. Let ^ := X]nm=o'^"."^ '4'n,m £ ^m(^^)- ^^ easy computation shows that 

+00 +00 

IkV'lliz = ^ \Cn,m\'^ n, ||aV||i2 = ^ \Cn,m\'^ (^ + 1) . 

n,m=0 n,m=0 

Since n^n + 1^2(n+|)^a(n+|) +^ holds true for any a > (arbitrarily small), 
then 

llaVlli. 



+00 


f 1^ 


2 +00 




a ^ Cn,ni ^ 


+ 2_^ \Cn,rn\ = 


= a\Hi^^\l2 +b\^/j\l2 


n,m=0 


\ ^y 


n,m=0 
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with b := ^ + \ where a" denotes either a or a^. I 

B Canonical transform for fast and slow variables 

This appendix is devoted to the concrete reahzation of the von Neumann unitary W intro- 



duced (in abstract way) in Section 4.2 The unitary W maps the fast and slow variables, 
which satisfy canonical commutation relation, into a set of canonical Schrodinger operators. 
In Section |B.l we derive a general version of the transform W "by hand", as a composition 



of three sequential transforms. In Section B.2 we compute the integral kernel of W 



r2J 



B.l The transform W built "by hand" 

Let l-i := L^(M^,(i^r) be the initial Hilbert space, with r := (r'i,r2). Let Q^ := {QniQ 
where Qr is the multiplication operator by Tj and Pj. := {Pr^.,Pr2) where P^ •= —ifid, 
with j = 1,2. Consider the fast and slow operators 

with a, /3 G C and f , w, v*,w* G M^ such that v ■ v* = w ■ w* = 1, v* ■ w = v ■ w* = and 
\v Aw\ = f > 0. 

Remark B.l. The choice v = b*, w = a*, a = JVq and /3 = JVq b defines the operators 




(56), while the choice v = v* = (0, —1), w = vf = (—1,0), a = /3 = 1 defines the kinetic 



momenta and the related conjugate operators introduced in Section |A.3[ ♦(> 



Observing that \a ■ Q,- ^ h ■ Pr]C ■ Qr ^ d ■ P^] = if^{a ■ d — b ■ cjl-^ one deduce that the 



operators (127) verify the following canonical commutation relations (CCR) 



[K^,K2]=^aHn, [G^,G2] = t/3Hn, [K,,G,]=0, i,j = l,2. (128) 

The Stone-von Neumann uniqueness theorem (see |BR97) Corollary 5.2.15) assures the ex- 
istence of a unitary map W {von Neumann unitary) 

W : n — >n^ := -Hs ® Hf := ^^(M, dx^) ® ^^(M, dxi) (129) 

such that 



WGiW-^ := Qs = multiplication by Xs, WGaW"^ := P, = -z/3^t^ (130) 



f) 
WKiW-^ := Q{ = multiplication by Xf, WiTsW"^ := Pf = -ia^^—. (131) 
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In other words, (Qs, -Ps) is a pair of operators which defines a Schrodinger representation 
on the Hilbert space "Hg := //^(IR, (ixg) while the pair (Qf, Pf) defines a Schrodinger represen- 
tation on the Hilbert space "Hf := //^(IR, dxi). Our purpose is to give an explicit construction 
for W. First of all consider the change of coordinates (ri,r2) i— )■ {ki := ^, k2 '■= ^). The 
inverse transforms are defined by ri{k) = |(u72^i — ^2^2) and r2{k) = j{vik2 — Wiki). The 
map J : L\R'^,£r) -^ L'^iB?,£k) defined by {J^lj){k) := ip{rik)) if ip & L^{R\dr) is 
unitary since the change of coordinates is invertible and isometric. Moreover JQr J~^ 
acts on l^^W? ^dfz) as the multiplication by rj[k), while JPt^J~^ = ^Pki + ^Pk2 where 
Pk '■= —iHdk, with j = 1,2. Then 

i (3'^ 1 i (3^ 1 

JGiJ := - Qfci —- Pk2, JG2J ■■= - Qk2 + -T-lj -Pfei (132) 

1 a i aB 1 1 a £ a/3 1 

JK,J-^ := --- Q,, - -f- P,„ JK2J-' := -- Qk, - -^j Pk,- (133) 

Let J^2,n '■ L^{R,dk2) — t- ^^(IR, ^^2) be the k2-Fourier transfor of weight fi, defined by 

(J^2,/.V')(C2) := \/M L^~'^''^''^'^(^2) d^^ and let Hi : L\R,dki) -^ L^{R,dCi) be the ki- 
parity operator defined by (112'?/') (Ci) •= "^(^Ci) (namely by the change of coordinates ki H- 
Ci). Let I the unitary map which identifies L'^{R'^,d'^k) with L^{R,dki) O L'^{R,dk2). Let 
Qc '■= {Q(ijQc2) whit Q^. the multiplication operator by Q and P^ := (P^j,P^2) where 
P^^ := —ifidc_i with j = 1, 2. One can check that 



Fix /i := — 2^2, then the unitary map £ := (Hi Cg) J-2,^ 



o . 



acts on the operators (127) in the following way 



-g^n niP,,nr^ = 


-Pa- 


i -^ L2(R 


,dCi)®L\\ 


(134) 

R,dC2) 




-Pc.) 


(135) 




+ PcJ. 


(136) 



CGiC '■— —^{Qci — Qci) ^ CG2J ■— — 

a / 
CK,C-':=--{Q^,+Q^,), CK2C'':= 

Now we can consider the change of coordinates (Ci,C2) '~^ {xs,X{) defined by 

^s = -^(Cl-C2) U' = --i{^^-^^^ 

^f= ^^(Cl + C2) [C2= ][^s+^X, 

The jacobian of this transformation is \d{Ci,C2)/d{xs,Xf)\ = ^1^1 =: G, then the map 
{TZilj){xs,Xf) := y/C il!{({xs,xi)) defines a unitary map TZ : ^^(IR^,^^^) — )■ L'^{R'^,dxs dx{). 
With a direct computation one can check that TZQc^TZ'^ acts on L'^{R?,dxs dxf) as the 

multiplication by Cj{xs,Xi), while nP(^^n~^ = ^ i^^^P^ + pX with j = 1,2. This 
shows that the unitary map W := I o 7^ o I~^ o C := I o TZ o 1^^ o (Hi (g> J^2,fi) o I o ^J is the 



von Neumann unitary that verifies the relations (130) and (131) 
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B.2 The integral kernel of W 

The unitary operators J : L'^{M?,(Pk) -^ L'^{R'^,d\) and 7^ : L'^{R'^,£C) ^ L'^iR'^,dxJxf) 
related to the change of coordinates (ri, r2) ^-> {ki, k2) and ((^i, (^2) '— ^ (^^sj ^^f) can be written 
an integral operators 



{Ji,){k)= J{r-k)i,{r)d\ {nv){x,,Xf)= R{C,x,,Xf)v{C)d\ 

with distributional integral kernels 

J{ri,r2;ki,k2) := 6 { n - -(^2^1 - V2k2) ] 6\r2 + -^{wiki - Vik2) 



R{CiX2\x,,x,) := VC 5 Ui + ji 



Xs xA\ 5 K2 - ^ ( a^s + - Xf 



with C" = ^ I f I • The /ci-parity operator Hi can be written as an integral operator with the 
distributional kernel 5{ri + ki) while the integral kernel of the A;2-Fourier transform J^2,ijli 

with/i:= -^, is 2^^*^^"'''. Then the unitary map I'^ o (Hi ®J'2,m) ol : L'^i^.d^k) -^ 
L^(]R^, (i^(^) as the integral distributional kernel 

L{k^MXxX2) :=5(Ci + A;i) -^^e*^«^'^ 

Summarizing the total transform W : L^(]R^,(i^r) — )■ L^(]R^,(iXs dx^) (up to the obvious 
identification I) can be expressed as an integral operator 



(WV^)(xs,Xf) = / W{r\x^,Xi) ^(r) d^r 
with a (total) integral distributional kernel 



W(ri,r2; Xs,X{) := — 1^=^= 6 { v ■ r — { Xs X{] ] e 2^^ ' " 

^/27r\a/3\ \ \ a ' ' 



(^ ^^ i^(^s+^ xA 
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